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Abstract 


In this thesis feasibility of the precise measurement of the Higgs-boson production 
cross section at the Photon Collider at TESLA is studied in detail. For the Standard- 
Model Higgs-boson production the decay to bb pairs is considered for the mass between 120 
and 160 GeV. The same decay channel is also studied for production of the heavy neutral 
Higgs bosons in MSSM, for masses 200-350 GeV. For the hrst time in this type of analysis 
all relevant experimental and theoretical effects, which could affect the measurement, are 
taken into account. 

The study is based on the realistic 77 -luminosity spectra simulation. The heavy quark 
background 77 — QQ{g) is estimated using the dedicated code based on NLO QGD 
calculations. Other background processes, which were neglected in the earlier analyses, 
are also studied: 77 — W^W~, 77 —r+r“, and light-quark pair production 77 — qq. 
Also the contribution from the so-called overlaying events, 77 — hadrons, is taken into 
account; a dedicated package called Orlop has been prepared for this task. The non¬ 
zero beam crossing angle and the hnite size of colliding bunches are included in the event 
generation. The analysis is based on the full detector simulation with realistic 6 -tagging, 
and the criteria of event selection are optimized separately for each considered Higgs-boson 
mass. 

In spite of the signihcant background contribution and deterioration of the invariant 
mass resolution due to overlaying events, precise measurement of the Higgs-boson pro¬ 
duction cross section is still possible. For the Standard-Model Higgs boson with mass of 
120 to 160 GeV the partial width r(h — 77 )BR(h ^ bb) can be measured with a statis¬ 
tical accuracy of 2.1-7.7% after one year of the Photon Gollider running. The systematic 
uncertainties of the measurement are estimated to be of the order of 2%. For MSSM 
Higgs bosons A and H, for M 4 = 200-350 GeV and tan/? = 7, the statistical precision 
of the cross-section measurement is estimated to be 8-34%, for four considered MSSM 
parameters sets. As heavy neutral Higgs bosons in this scenario may not be discovered 
at LHG or at the hrst stage of the collider, an opportunity of being a discovery 

machine is also studied for the Photon Collider. 
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Chapter 1 
Introduction 


So far, all experimental results concerning fundamental particles and their interactions 
are well described by the Standard Model (SM), consisting of the electroweak theory 
(EWT) and the quantum chromodynamics (QCD). These theories allow us to quantify 
electromagnetic, weak and strong processes; among all known kinds of interactions only 
the gravity is not incorporated in the SM framework. The very important ingredient 
of the SM is the so-called Higgs boson, h, which is responsible for generating masses 
of all particles. Although predicted by the model, the Higgs boson has not yet been 
experimentally detected. However, such a particle must exist if the SM is to remain a 
consistent theory. Consequently, a search for the Higgs boson is among the most important 
tasks of the present and future colliders. Once the Higgs boson is discovered, it will be 
crucial to determine its properties with high accuracy, to understand the mechanism of 
the so-called electroweak symmetry breaking (EWSB). 

The neutral Higgs boson couples to the photon pair only at the loop level, through 
loops of all massive charged particles. In the SM the dominant contribution is due to W 
and t loops. This loop-induced /177 coupling is sensitive to contributions of new particles 
which may appear in various extensions of the SM. Hence, the precise measurement of the 
Higgs-boson partial width T{h ^ 77 ) can indicate existence of very heavy particles even 
if their direct production is not possible. A photon-collider option^ of the collider 

offers a unique possibility to produce the Higgs boson as an s-channel resonance in the 
process 77 — h. As the SM Higgs boson with the mass^ below ~ 140 GeV is expected 
to decay predominantly into the bb final state, we consider the measurement of the cross 
section for the process 77 ^ h —> bb, shown in Eig. II .11 for the Higgs-boson mass in 

photon collider option was foreseen for all projects of the e+e" linear collider: TESLA Q], NLC 
[2 and GLC (earlier JLC) OBI- In this work the photon collider at the TESLA is considered. The 
superconducting technology developed within the TESLA project has been recently selected as the best 
suited for the International Linear Collider. Decision of the International Technology Recommendation 
Panel was presented during the ICIIEP2004 conference in Beijing 0. 

^The energy unit [GeV] is used for masses and momenta, i.e. the speed of light is set to 1. However, 
for lengths and times the corresponding units are [m] and [sj. 
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Figure 1.1: A diagram of the process 77 —/r — bb. The Higgs boson particle couples to 
photons through the loop of all massive and charged particles. 


the range = 120-160 GeV. The aim of this study is to estimate the precision with 
which this measurement and extraction of r(/i — > 77 ) will be possible after one year of 
the TESLA Photon Collider running. 

Besides precision measurements, a photon collider can be also considered as a candi¬ 
date for a discovery machine. In case of the Minimal Supersymmetric extension of the 
SM (MSSM) the photon collider will be able to measure the production cross section of 
the heavy neutral Higgs bosons, A and FT, covering the so-called “LHC wedge” in the 
MSSM parameter space, i.e. region of intermediate values of tan/?, tan/? 4-10, and 
masses Ma,h above 200 GeV. For this part of parameter space, MSSM Higgs bosons A 
and H may not be discovered at the LHC pmiH) and at the first stage of the e+e 
linear collider [Hj becanse of small branching ratios into leptons or photons (which allow 
the efficient signal selection) and because of the kinematical limit Ma,h ^ for 

pair production process AH, respectively. Parameter range considered in this 

analysis corresponds to a SM-like scenario where the lightest MSSM Higgs boson h has 
properties similar to the SM Higgs boson, while heavy neutral Higgs bosons are nearly 
degenerated in mass and have negligible couplings to the gauge bosons W/Z. We consider 
the process ■yj ^ A, H —>■ bb at the Photon Collider at TESLA for Higgs-boson masses 
Ma = 200-350 GeV. The aim of the presented study is to evaluate the discovery poten¬ 
tial of the considered experiment by estimating the statistical significance of the signal 
measnrement for the chosen region in the MSSM parameters space. Also the precision of 
the yy ^ A, H ^ bb cross section measurement is estimated. 

The measurements of T{h —> 77 )BR(fi —> bb) and a (77 A, H ^ bb) at a photon 
collider have already been stndied before [iniiniiiaiiHiiiaiiaiiniiniiiHiiiniiaEniini 
EaEHlElEniEniEIll^HlEllIHQllHI]. Although very promising estimates were obtained, 
many important aspects of the measurement were not considered. This study is the first 
one to take all relevant experimental and theoretical effects into account. Only results of 
such a realistic analysis can be used to support the project of the Photon Collider in the 
framework of the International Linear Collider. 

The motivation for this study is outlined in Chapter [21 The proposed experimental 
setup and simulation tools are described in Chapter |H1 In Chapter |31 details of the signal 
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and background simulations are given. A discussion of the event selection and the final 
results for SM and MSSM scenarios are given in Chapters El and IHl respectively. All results 
presented in Chapters El El and El and in Appendices were obtained by the author of this 
thesis. 
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Chapter 2 
Motivation 


In this chapter our current understanding of the Higgs mechanism and prospects for a 
discovery of the Higgs boson are outlined. The Higgs sector in the SM is discussed first, 
then its extension to the MSSM is shortly reviewed. Current limits on the Higgs-boson 
mass from direct and indirect measurements are summarized. Expected experimental 
results at future colliders, relevant to the presented study, are also given. An in-depth 
review of the Higgs-boson theory and phenomenology can be found, for example, in |32[ 
EB]- An extensive summary of experimental results concerning Higgs-boson searches is 
presented in Q. 


2.1 The Higgs sector in the SM 

Among all fundamental particles of the SM only the Higgs boson still remains hypothet¬ 
ical. This neutral spinless particle is required in the model to break the gauge symmetry 
of weak interactions. Photon, which is a carrier of electromagnetic force, is massless. But 
three weak bosons Z, kT+ and W~ are massive; this is a serious difficulty as SM equa¬ 
tions for interactions involving massive bosons lack a very basic property, the so-called 
gauge invariance^. Other problem emerges in cross section calculations for some weak 
processes, e.g. e+e“ —> W^W~, because unitarity condition is violated for this transi¬ 
tion. Probability current is not conserved unless we introduce new particles which couple 
to electrons and massive bosons. One complex Higgs doublet (four real scalar fields) 
is introduced in the SM in order to describe experimental results and to preserve clear 
theoretical picture. These new fields, filling the vacuum, couple to the massless vector 
bosons, giving them effective mass. This mechanism, introduced by P. Higgs ra, allows 

^ A principle of gauge invariance originates in the classical theory of electromagnetism and reads: there 
is a transformation of electromagnetic four-potential after which physically relevant fields, E and B 
(or the tensor = d^A'^ — d’'A^), remain unchanged. In quantum field theories the invariance of 
equations after simultaneous, special transformations of all fields is required. Thus, the gauge invariance 
principle determines the allowed interaction terms. 
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US to introduce massive gauge bosons in the theoretical description without violating the 
gauge invariance (so-called spontaneous symmetry breaking). One of the scalar fields is 
expected to exist as a real particle, so-called Higgs boson, h. All conplings of the Higgs 
boson to other particles and its self-couplings are predicted by the SM; the couplings to 
bosons (fermions) are proportional to the mass squared of the boson (the mass of the 
fermion). The only nnknown parameter of the theory is the Higgs-boson mass, M^. An 
intelligible introdnction to the Higgs mechanism can be fonnd, for example, in jdfij . 

The SM constitntes a complete effective theory of fnndamental interactions (exclnding 
gravity). Existence of the new particle, the Higgs boson, explains how the electroweak 
symmetry (or gange invariance) is broken and solves the nnitarity problem in weak re¬ 
actions. However, this great theoretical achievement is nndermined by some nnsolved 
problems. On the way to the Planck energy scale some new phenomena are expected to 
appear. Otherwise, without unnatural tuning, higher order corrections to the Higgs-boson 
mass diverge as the energy scale increases (so called “hierarchy problem”). The second 
problem is due to our expectation that at some high energy scale all interactions should 
unify («.e. their couplings should be equal) which is not exactly the case in the SM.^ To 
fulfill this unification reqnirement new particles or interactions have to be introdnced. 


2.2 Higgs sector in the MSSM 

The new symmetry between bosons and fermions, so-called supersymmetry (SUSY), conld 
remove the two above-mentioned problems of the SM. It gnaranties cancellation of diver¬ 
gences in Higgs-mass calculation.^ Also the unification of three fundamental couplings 
is realized. However, in the general case of the Minimal Supersymmetric extension of 
the SM (MSSM) aronnd 100 new parameters mnst be introdnced whose valnes are not 
predicted by the model. All SM particles have their superpartners: fermions - spin-zero 
bosons {e.g. electron - selectron), bosons - fermions (e.g. higgs - higgsino, W - lU-ino, 
photon - photino). To generate masses for all particles and sparticles, two Higgs donblets 
{i.e. eight fields) have to be introdnced. As a result, supersymmetric models contain five 
Higgs bosons (instead of one Higgs particle). Two of them are neutral scalars and are 
denoted as h and There is also one nentral psendoscalar. A, and two charged scalars: 
H+ and R-. 

The Higgs sector of the MSSM is described by a subset of parameters which includes: 

^Only approximate unification is obtained in the SM. At the scale of Mqut ~ 10^® GeV couplings are 
’unified’ to 0(10%) UTj . 

^In fact, cancellation is not exact as the supersymmetry is broken, i.e. particles have different masses 
than their superpartners. This results in the prediction that masses of superpartners cannot be heavier 
than a few TeV. Otherwise supersymmetry does not solve the hierarchy problem. 

"^By definition, h denotes the lighter scalar Higgs boson, and H denotes the heavier one. 
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1. tan/3 - the ratio of vacuum expectation values of neutral Higgs fields coupling to 
up- and down-type fermions, tan/3 = Vu/vd- 


2. Ma - the mass of the neutral, pseudoscalar Higgs boson, A. 

3. /(i - the supersymmetry-breaking higgs-higgsino mass parameter. 

4. M 2 - the supersymmetry-breaking universal gaugino mass parameter (mass of the 
kb-ino; masses of other gauginos are related with M 2 ). 

5. Mj^, Mj^ and Aj - other supersymmetry-breaking parameters: masses of left- and 
right-handed supersymmetric partners of fermion / and their coupling to Higgs 
bosons, respectively. 

Only first two parameters, tan/3 and Ma, infiuence the Higgs sector on the tree level. 
Other parameters can affect properties of the Higgs bosons via radiative corrections. In 
contrast to the SM, mass of the lightest Higgs boson, h, is constrained, i.e. h cannot be 
heavier than around 150 GeV. 


2.3 Status of the Higgs-boson searches 

In precise calculations of the SM predictions the higher order corrections resulting from the 
Higgs boson contribution are sizable and must be taken into account. Expected results 
for many observables depend on the Higgs-boson mass, M^. Thus, constraints on the 
value of Mh can be obtained from the analysis of electroweak measurements. The result 
of such analysis is shown in Fig. 12.11 pkH] . where the value from the SM fit to precise 
measurements at LEP, SEC, Tevatron and other experiments is presented as a function 
of the Higgs-boson mass. The best agreement is found for M^ ^ 126 GeV, and with 95% 
G.L. the upper limit on Mh is 280 GeV. The best fit value is slightly above the lower mass 
limit from the direct searches at LEP; excluded is the mass range Mh < 114.4 GeV [221 • 
In the MSSM case limits for the Higgs-boson masses depend on other model param¬ 
eters. In the general approach, when other parameters are allowed to vary, we can only 
conclude that all Higgs bosons must be heavier than 80-90 GeV if model with no CP- 
violation in Higgs sector is assumed. However, some MSSM parameter sets result in the 
lightest higgs, h, having couplings similar to those of the SM Higgs boson (SM-like sce¬ 
narios). In such cases the Higgs-boson mass constraints are similar to those obtained in 
the SM. 
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Figure 2.1: The result of the SM fit to the precision electroweak measurements jnHj- The 
~Xmin value is presented as a function of the Higgs-boson mass (here denoted as 
mn)- The best agreement is found for Mh ~ 126 GeV, and with 95% C.L. the upper limit 
on Mh is 280 GeV. Excluded mass range from direct searches at LEP, Mh < 114.4 GeV, 
is indicated in yellow. The band represents the theoretical uncertainty due to missing 
higher order corrections. 


2.4 Prospects for Higgs-boson measurements 


If the mass of the SM Higgs boson is around 115 GeV it is still possible that it will 
be discovered at the Tevatron. However, only future machines will have sufficient higgs 
production rates to measure precisely the mass and couplings of the Higgs boson(s). All 
large accelerator projects aim at measurements of the Higgs-boson properties from which 
the fundamental one is the mass, Mh, being at the same time the only unknown parameter 
in the SM. Measurements of other parameters describing the Higgs boson (total and partial 
widths, branching ratios, spin, parity) are considered as the crucial tests of the SM and 
its extensions. Below the estimated precisions of future measurements are summarized, 
in the expected order of appearance. 
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Large Hadron Collider 

At the Large Hadron Collider (LHC), which should become operational in 2007, Higgs 
boson(s) will be produced in processes of gluon fusion (about 80% of the SM Higgs-boson 
production rate for Mh < 2Mz) and W boson fusion. Both experiments, ATLAS and 
CMS, have presented detailed studies showing that the SM Higgs boson will be discovered 
at LHC if it is lighter than about 1 TeV. The Higgs-boson mass can be measured with 
precision about 0 . 1 % for Mh < 400 GeV [HI] . 

Various ratios of Higgs partial widths can be determined with precisions about 10-20%, 
assuming integrated luminosity of 100 fb“^ jH]. For heavy SM higgs, Mh > 250 GeV, 
also its total width, F/j, can be measured as it becomes larger than the experimental mass 
resolution n. The heavy MSSM Higgs bosons, A and FT, will be observed at the LHC 
for most of the allowed MSSM parameter space. However, there is a region of M 4 — tan/I 
values where the LHC may not be able to discover heavy MSSM Higgs bosons. This is 
the so called “LHC wedge” covering Ma > 200 GeV and tan/3 4-10. 

Linear e+e“ Collider 

According to the currently proposed schedule, the International Linear Collider 

(ILC) can become operational in 2015. Two processes contribute to the Higgs-boson 
production at the ILC: Higgs-strahlung and vector boson fusion. For the SM Higgs 
boson with mass in the range Mh ~ 115-180 GeV the expected precision of the mass 
measurement can be better than 0.05% jH- As the background is much smaller than at 
the LHC, Higgs-boson branching ratios can be determined with much better precision 
and in the model-independent way. Branching ratios BR(h — 77 ) and BR(h — bb) may 
be determined with accuracy of about 10% and 1.5%, respectively, after one year of ILC 
running at nominal luminosity EaESl. 


Experiments after LHC and LC 

After LHC and ILC measurements there will still be some properties of the Higgs-boson(s) 
which are poorly known and should be determined with greater precision at other experi¬ 
ments. One of the interesting quantities is the 77 partial width, F(h — 77 ) = F/^ BR(h — 
77 ). As already mentioned in Chapter [U measurement of partial width F(h —> 77 ) is of 
special importance as it is sensitive to all charged particles which have mass generated 
by the Higgs mechanism. Due to the so-called non-decoupling effect contributions of the 
heavy charged particles to the h —> 77 loop are hnite even in the limit of inhnite mass 
of the particle. Thus, the measurement of F(h — 77 ) can indicate existence of particles 
whose direct production will be impossible in available accelerators. 
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Figure 2.2: Ratio of the higgs —*• 77 decay widths in the 2HDM and the SM as a 
function of the higgs mass, Mh,H, for solutions A (left) and B (right) jl2|. The yellow 
bands correspond to the uncertainties expected at Mh,H = 120 GeV. 


In the MSSM, the loop indnced coupling is sensitive to contribntions of snpersymmetric 
particles EH EEl ESI: chargino and top sqnark loops can lead even to 60% difference 
between the SM and the SUSY couplings. Scenarios, in which all new particles are very 
heavy, may be realised not only in the MSSM but also in other models with extended 
Higgs sector, for example in the Two Higgs Donblet Model (2HDM). In this case the two- 
photon width of the Higgs boson will differ from the SM value due to the contributions of 
the heavy charged Higgs bosons, even if all direct couplings to gauge bosons and fermions 
are eqnal to the corresponding SM conplings. Different realizations of the 2HDM have 
been discnssed in B3- Assuming that the partial widths of the observed Higgs boson to 
quarks, Z or W bosons are close to their SM values (SM-like scenario), three different 
combinations of conplings are possible. Fig. l2.2l shows deviations of the two-photon Higgs 
width from the SM valne for the three SM-like solntions considered in na. For solntion 
Bu one expects significant deviation from the SM predictions since, as compared to the 
SM, there is a change of the relative sign of the top-quark and the W contributions. 
Consequently, for solntion Bu the higgs —> 77 width is significantly larger than in the 
SM, where these two contributions partly cancel each other. However, deviation due to 
the charged Higgs-boson contribution only (solution A) is much smaller, of the order of 
5-10%, and the measurement with precision at the level of a few percent is reqnired. 

The machine best suited for the measurement of the Higgs boson two-photon width 
is a photon collider. All photon collider proposals (within TESLA, NLC, GLC and CLIC 
projects) emphasize the feasibility of a very precise F(h — 77 ) measnrement. For the 
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light SM-like Higgs boson the most promising process is 77 —/r —66 due to the very 
high branching ratio BR( 6 , —> 66 ). The measurement of the cross section for this process 
has been studied in detail and is the main subject of this work. One has to note that the 
final results on r(fi — 77) from a photon collider must rely on BR(/r — 66 ) measurement 
from other experiment - this branching ratio should be determined at LC with precision 
of around 1.5% IlH]. Measurements of r(6, —> 77)BR(6, —> 66 ) at the Photon Collider and 
of the branching ratios BR(h —77) and BR( 6 , — 66 ) at the collider can be used to 
determine the total Higgs width, P/^, in a model independent way. 

The quantity r(6, ^ 77 )BR(fi — 77 ) could also be unfolded from the cross section 
measurement for the process 77 ^ 6 . —> 77 . Unfortunately this process has a very low rate 
due to small 77 branching ratio. Moreover, in addition to a one-loop background process 
77 77 also ’machine’ background 07 ^ 07 must be considered. Even with optimistic 

assumptions about angular coverage (down to 3°) and high granularity of calorimeter the 
precision of the measurement has been estimated to be of the order of 30% uni- Therefore 
this measurement cannot be considered as an alternative to the analysis of the process 
77 —> 6, —> 66. 

The Photon Collider (PC) seems to be the only machine allowing precise measurement 
of r(/r — 77). However, physics potential of the PC is much reacher and complementary 
to that of the LHC and the ILC. For Higgs sector itself many interesting measurements 
can be considered: 

• CF-parity of the Higgs boson. The CP-parity of the Higgs boson can be determined 
in a model independent way from analysis of angular correlations in 4-fermion decays 
h — W^W~ —4j and h ^ ZZ ^ 2j 21 [IHl- Similar measurements have been also 
proposed for LHC and ILC. 

• Phase of the 6, —> 77 coupling. In addition to T{h —> 77), precise measurements 
at the PC are also sensitive to the phase of the 6, —> 77 amplitude, (j)h-y-y The 
phase (j)h-y'y can be extracted from the measurement of the interference between the 
resonant higgs production processes 77 — 6, — W^W~ and the background process 
77 —> W^W~ |1D]. It turns out that for higgs mass of the order of 300-400 GeV the 
phase (ph-y'y is more sensitive to possible contributions of heavy charged scalar particle 
than r(6, —77). Only by combining r(6, — 77) and (ph'yy measurements at the 
PC with those at the LHC and at the LC unique determination of the Higgs-boson 
couplings and distinction between various models will be possible. 

• Charged Higgs-bosons production. For intermediate values of tan f3 the charged MSSM 
Higgs boson, , may not be discovered at the LHC, if its mass is greater than 
~ 150 GeV [3 IH]- In the photon collider H^H~ pairs can be produced in QED 
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process and the mass reach for discovery can be extended up to ~ 300 GeV (if run¬ 
ning at -y/s^ = 800 GeV). At the LG, running with = 800 GeV, almost 
an order of magnitude smaller number of higgs pairs could be produced during the 
same time due to smaller cross section pm. 

The second generation linear collider project for which the feasibility study is 

still in progress is CLIC ra. At this accelerator multi-TeV energies will be obtained, and 
the wider range of masses will be accessible for new-particle searches. The photon-photon 
collision mode has been proposed for this machine as well, but the detailed physics case 
studies are still missing. 

Opportunity of studying elementary particle collisions at midti-TeV energies is the 
main reason for considering the next generation project of a muon collider. As muons 
lose much less energy via bremsstrahlung than electrons a circular accelerator option is 
preferable even for much higher beam energies than those accessible at LEP2. Circulating 
beams would allow the operation with high luminosity. However, progress must still be 
made in formation of high-intensity /x-beams with low emittance. At the muon collider 
the same channels could be used for Higgs-boson measurements as at collider but 
with extended mass reach. Moreover, the s-channel production in the process p/x ^ higgs 
has cross section sufficient for precision measurements due to higher mass of the muon 
0 . As the energy spread of muon beam is expected to be negligible, energy scan at 
the muon collider would result in determination of the Higgs-boson mass and width with 
precision of the order of 2 MeV. 

Construction of a photon collider based on the muon collider is also possible. Unfor¬ 
tunately the high mass of the muon is a problem in this case. The maximal energy of 
photons from Compton backscattering would be very small compared to the beam en¬ 
ergy, e.g. the maximal photon energy would be only around 75 GeV at 20 TeV /x-beam as 
compared to about 400 GeV for 500 GeV electrons with the same laser setup^. 


®Laser parameters of the TESLA Photon Collider design are assumed. 
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Chapter 3 


Collider and detector 


3.1 The TESLA Linear e+e Collider 

Results presented in this thesis are based on the design and machine parameters of the 
TESLA (TeV-Energy Superconducting Linear Accelerator) collider [221 • Accelera¬ 

tor design is based on the superconducting technology, recently accepted by International 
Technology Recommendation Panel as the best suited for the ILC project |2]. Each of 
the two linear accelerators, which accelerate e~^ and e~ towards the interaction region, 
will consist of around 10000 one-meter-long superconducting cavities. Cavities made from 
niobium and cooled to temperature of 2 K can provide accelerating field with gradient well 
above 35 MV/m. As the average gradient equal to 23.4 MV/m is required for operation 
at the nominal total collision energy of 500 GeV, the opportunity emerges to increase the 
machine reach up to 800 GeV or even 1 TeV. With dumping rings and other additional ac¬ 
celerator components the total length of the machine is 33 km. Main machine parameters 
are listed in the table im 


Description of the parameter 

Value of the parameter 

Accelerating gradient 

No. of accelerating structures 
Train repetition rate 

No. of bunches per train 
Bunch spacing 

No. of e“(e’'') per bunch 

Beam size at IP (cr^iCTy) 

Bunch length at IP (cr^) 
Luminosity 

Luminosity per year 

23.4 MV/m 

21024 

5 Hz 

2820 

337 ns 

2 ■ 10^° 

553 nm; 5 nm 

0.3 mm 

3.4 • 10^^ cm“^s“^ = 34 nb“^s“^ 
340 fb- 


Table 3.1: Main parameters of TESLA Linear Collider for the 500 GeV design. One 
accelerator year [y] is assumed to be equal to 10^ s. 
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3.2 A photon collider as an extension of the LC 


Future linear colliders offer unique opportunities to study photon-photon inter¬ 

actions if the idea of a photon collider is realized. In this option the energy of the 
primary electron-electron^ beams is “transfered” to photons in the process of Compton 
back-scattering Assuming the beam electron collides with one laser photon head-on, 
the highest energy of the scattered photon, ^ is eqnal to: 

pmaxl _ ^ E^ + Pe 

'y 9 ’ 


where 


X = 


2EL{Ee+Pe) 


mi 


El, Ee, pe and are the energy of laser photon, the energy of beam electron, its 
momentum and mass, respectively. As the energy of primary electrons will be of the 
order of 100 GeV, one can use an approximation pe ^ E^. In this case the simplified 
formnlae are obtained: 


E, 


maxi 

7 


X 


X = 


X +1 

^EiEe 


-E, 


mt 


Collision of a high energy photon from Compton back-scattering with a laser photon in 
the conversion region can result in creation of pairs. This process would significantly 
limit the 77 -lnminosity as its cross section is comparable with that of the Compton 
scattering. Therefore, it is preferable to select laser parameters snch that the threshold 
for pair creation is not reached. The condition which has to be imposed on the 

invariant mass of two photons is: 


< 2m,. 

This is equivalent to the requirement: X < Xthr = 2 + 2\p2 K. 4.83. 

Leading order Compton cross section formula indicates that the approximate monochro¬ 
maticity of the 7 beams can be achieved for high polarizations of colliding electrons and of 
laser photons. Figures ini and lT^ show the expected distribntions of prodnction probabil¬ 
ity and of the circular polarization of high-energy photons, as a function of the photon 
energy relative to the primary electron energy, y = E^y/E,, for various combinations of 
laser photon polarization, P,, and electron polarization, P,. With the choice P, = 1 and 
Pc = — 1 most of the back-scattered photons have energies close to and are highly 

^For the photon collider positron beam is not needed. Moreover, use of two electron beams has 
important advantages: higher photon polarization and reduced backgrounds in the interaction region. 
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Figure 3.1: Probability distributions for production of photon with relative energy 

y = E^/Eg. Results for various combinations of laser and electron polarizations, Pc and 
Pe, are shown. 

polarized. 

In the laboratory frame photons from Compton back-scattering are strongly boosted 
in the electron-beam direction. The angular distribution has the width corresponding to 
the characteristic angle 9ch ~ l/ 7 e where je is the relativistic factor for the beam electron, 
7 e = Eel mg. For Ee = 250 GeV one obtains 9ch ~ 2 /rrad, i.e. the photon beam is strongly 
collimated along the incident electron beam direction. Therefore, two high energy photon 
beams produced in Compton back-scattering can be collided head on in the setup shown 
in Fig. 13.31 As the distance between conversion point (CP) and interaction point (IP), 5, 
will be of the order of 3 mm, one can see that the additional spread of the photon bunch 
will be only about 6 nm which is much smaller than the transverse size of the e“-bunch in 
the X direction and comparable with the e“-bunch size in the y direction. Therefore, the 
77 -luminosity will be of the same order of magnitude as the e“e“ geometrical luminosity. 

Studies of the effects present in the conversion and interaction points revealed that 
also the following corrections should be included in a description of 77 -luminosity spectra: 

1 . Correlations between energy and scattering angle of Compton photons. As more 
energetic photons scatter with smaller angles, high-energy photons in the ’core’ 
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Figure 3.2: Circular polarization of scattered high energy photons, as a function of 
relative photon energy, y = E^/Ee- Results for various combinations of laser and electron 
polarizations, Pc and Pg, are shown. 

of the beam collide with high-energy photons of the opposite beam with greater 
probability than with low-energy photons forming beam ’halo’. Thanks to this effect 
high-kF.),.^ part of 77 -luminosity is enhanced in comparison to the simple convolution 
of both spectra [HHl- 

2 . An effective increase of the electron mass due to its transverse motion in the strong 
electromagnetic field of the very intense laser beam: ml —> rUe// = ml{l -|- 
Here the ^ parameter is related to the strength of the electromagnetic field in the 
conversion region and is used to describe nonlinear effects. 

3. Scattering of electrons on two laser photons: e~ -|- 27 /, —^ e“ -|- 7 . Interactions with 
three and more laser photons are snpposed to be negligible. 

4. Interactions of laser photons with electrons which already scattered one or more 
times. 

5. Nonlinear pair creation 7 -f myi e'^e~ which should be taken into account 
even for x < Xthr- 
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Figure 3.3: A basic scheme of the photon collider (HU. Primary electrons scatter on the 
laser photons (region C) at the distance b from the interaction point (IP). High-energy 
photons prodnced in the Compton back-scattering follow the electron-beam direction and 
collide with photons coming from the opposite side. 


6 . Coherent e+e pair creation by a high energy photon in the electromagnetic field 
which is present in the interaction point. 


Before the interaction point, to minimize some of aforementioned effects, the possibility 
was studied to remove electrons from the beam with special magnets. For a design with 
6 ~ 1 cm a small magnet was foreseen with magnetic field B = 5 kC deflecting electrons 
before the IP. However, in the current design, optimized for highest luminosity, this is no 
longer possible dne to the short distance of 3 mm between conversion point and IP. 

The finite beams-crossing angle at the interaction point, with “crab-wise” tilted elec¬ 
tron bnnches PI has been recently accepted as the solntion for the linear collider. This 
is a preferred scheme for a photon collider because the removal of high-energy-photon 
bunches after the interaction would be very difficult with collinear beams. Crab-crossing 
solntion preserves the same luminosity as for head-on collisions. However, electromag¬ 
netic interaction between beams must be included in the full simulation of 77 -luminosity 
because primary electrons are traversing through interaction point. 

The more complete description of processes outlined here and other effects influencing 
77 -lnminosity spectra can be fonnd in |55] . 
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3.3 The Photon Collider at TESLA 


According to the current design of the Photon Collider at TESLA jTj, the energy of 
the laser photons is assumed to be fixed for all electron-beam energies. Laser photons 
are assumed to have circular polarization Pc = —1, while longitudinal polarization of 
electrons is Pg = 0.85. This configuration of polarizations corresponds to the energy 
spectra of back-scattered photons peaked at high energy (see Fig. HH). With the same 
choice of parameters for each beam we maximize probability that two high-energy photons 
will collide with the same polarization, i.e. in the state with total angular momentum, J, 
equal to zero, so a spinless resonance can be produced. 

To profit from the peaked 77-luminosity spectra the energy of primary electrons has 
to be adjusted in order to enhance the resonance production signal at a particular mass. 
The use of a by-pass for electron beams is considered if the energy much lower than the 
nominal one is required. In this case the luminosity will be approximately proportional 
to the beam energy, L oc Pg- 

3.3.1 Photon-photon luminosity spectra 

As described in Ea the 77-luminosity spectrum is influenced by many various effects. 
To take them properly into account, a dedicated program for detailed beam simulation 
for the Photon Collider at TESLA has been developed isni. Large samples of 77 events 
were generated at selected energies and are available for further analysis. The simulated 
photon-photon events were directly used in this analysis when a proper description the 
low energy tail of the spectrum was crucial, e.g. for the so-called overlaying events'^. 
However, in the high-energy part of the 77 spectrum, i.e. for lTy.y > 0.3 where 

the results of the full simulation are well described by the CompAZ 
parametrization The subroutines of the CompAZ package were used when “con¬ 

tinuous” description was necessary. For example, in case of a very narrow resonance 
production the full simulation provides only a few 77 events in the region of interest per 
one million of simulated photon-photon collisions. Hence, analytical approach is much 
more efficient. Also the NLO QCD program used for generating 77 — QQ{g) events 
required a functional description of the luminosity spectrum for a proper calculation of 
the cross section. 

This analysis is based on the realistic 77 luminosity simulation for the Photon Collider 
at TESLA |Tj. Some earlier studies of Higgs-boson production in the process 77 —> 
higgs bb assumed other laser parameters and/or “ideal” 77-luminosity spectrum (i.e. 
spectrum corresponding to the LO Compton cross-section formula) [221 IIBI 121] ■ The 
“ideal” spectrum, used in is compared with CompAZ parametrization of the realistic 
^See Section in Chapter 01 
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Figure 3.4: Photon-photon luminosity spectra used in the analysis of the SM Higgs- 
boson production with mass Mh = 120 GeV, as a function of the invariant mass of the 
two colliding photons The spectrum obtained from CompAZ parametrization is 

compared with the lowest order QED predictions for the Compton scattering, used in the 
earlier analysis (lines). The total luminosity distribution (J^ = 0, ±2) and the 7^ = 0 
contribution are shown, separately. 


beam simulation results in Fig. 13.41 As can be seen, the “ideal” spectrum would be 
more advantageous for the narrow-resonance production. Additional effects, which have 
to be taken into account in the realistic study, increase the contribution of low energy 
77 collisions and make the high energy peak wider. Moreover, the leading order results 
for Compton process assume a; = 4.8 whereas fixed laser wave length is assumed in the 
present design for the whole energy range of electron beams, resulting in x parameter 
values smaller than the optimum value used in an “ideal” spectrum. Therefore, results 
obtained in this analysis, using the realistic spectra description, should not be directly 
compared to results obtained with “ideal” spectra. 

Details of the kFy..y-spectrum obtained with CompAZ for = 210 GeV are shown in 
Fig. 13.51 Contributions from two polarization combinations are indicated separately, i.e. 

= 0 and \Jz\ =2 where Jz is the total 77 angular momentum projected on a collision 
(z) axis.^ The suppression of = 2 luminosity can be clearly seen in the high-kFy.y 

^The coordinate system used in this document is a right handed system, with the z-axis pointing in 
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part of the spectrum. The threshold at = 131 GeV, expected for collisions of two 

photons produced in the lowest order Compton scattering, is not sharp. There is a tail 
of collisions involving two photons from the second order process e~ + 27 ^, —e“ + 7 for 
which the highest possible energy is around 161 GeV. The intermediate structure emerges 
from “mixed” collisions with photons originating from different scattering processes. 

If not stated explicitly otherwise, the resnlts presented in this work are obtained for an 
integrated luminosity expected after one year of the TESLA Photon Collider running inn]. 
In Table ESm the total photon-photon luminosity per year, is shown for different 
electron beam energies. Also shown are: the Higgs-boson mass corresponding to the max¬ 
imum of J 2 = 0 luminosity spectrum for given beam energy, and the expected luminosity 
in the high energy part of the spectrum, i.e. for > IT^™ where IT™“ = 0.5 . 


3.3.2 Collision region 

The e“e“ beams used in the photon collider will have similar geometrical parameters as the 
beams of LG (they will be prodnced in the same damping rings, compressed by the 
same compression system etc.). However, beamstrahlung due to beam-beam interactions 
is not present and the e“e“ beams can be focnsed on a smaller area at the interaction 
point (IP). For ^/s^e = 200 GeV electron bnnches are assnmed to have: = 140 nm, 

ay = 6.8 nm and az = 0.3 mm. The longitudinal (z) photon-beam bunch size at the 
photon collider is approximately the same as the corresponding size of an electron bunch. 

the direction of the electron beam in the e+e" mode, and the y-axis pointing upwards. The polar angle 9 
and the azimuthal angle 4> are defined with respect to z and x, respectively, while r is the distance from 
the z-axis. When describing selection procedure the angle with respect to the beam direction is limited 
to 0 < 0 < n/2. 


[GeV] 

Muggs [GeV] 

L.y.y [fb- 1 ] 

IT““ [GeV] 

L.,^(IT^^ > IT””) [fb-i] 

211 

120 

410 

65 

111 

222 

130 

427 

70 

116 

234 

140 

447 

75 

121 

247 

150 

468 

81 

126 

260 

160 

489 

86 

132 

305 

200 

570 

106 

150 

362 

250 

683 

131 

173 

419 

300 

808 

157 

196 

473 

350 

937 

182 

216 


Table 3.2: Luminosity per year expected for the Photon Collider running at given y/s^. 
One accelerator year [y] is assumed to be equal to 10^ s. Also shown are: the Higgs- 
boson mass, Mhiggs, corresponding to the maximnm of the = 0 Inminosity, and the 
expected Inminosity in the high energy part of the spectrnm, i.e. for ITy-j, > IT™“ where 

y^min ^ g_5p^maxl_ 
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Figure 3.5: CompAZ parametrization of the photon-photon luminosity spectra for 

= 210 GeV as a function of the invariant mass of two colliding photons The 

total luminosity and contributions with J^=0 and = ±2 are shown. 


However, as a distance of 6 = 2.6 mm is foreseen between CP and IP, transverse sizes of 
the photon bunch will be greater than that of the electron bunch due to the angular spread 
of the Compton scattering. This affects distribution in the y direction, as an additional 
spread is of the order of ay, but it does not influence the x-size of the bunch (ax 3> ay). 

For two head-on colliding bunches, which have Gaussian distribution with equal vari¬ 
ances (and the same speed), the spacial distribution of collision probability follows 3- 
dimensional Gaussian distribution with all three variances two times smaller than corre¬ 
sponding bunch parameters, i.e. a'^ = a\l2 (a = x,y,z). 

Transverse dimensions of a photon bunch decrease slightly with For ^/^ = 

200-800 GeV x and y dispersions of the photon bunch are 140-70 nm and 15-5 nm, 
respectively. Vertical dimension of IP density, a'y = ay/\/2, is about 10 nm or smaller, 
so distribution in this direction is too narrow to influence the event reconstruction and 
can be safely neglected. So would be the horizontal dimension a'^ if the beams collided 
head-on. However, the crab crossing scheme results in modified collision density in the 
horizontal direction. This effect is schematically shown in Fig. 13.61 Assuming that beams 
collide with relative angle Oc = 34 mrad, the x-size of collision region is given by the 
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Figure 3.6: A scheme of the crab-wise bunch-bunch collision with the relative angle Oc 
between beams. Bunches (the blue ellipses) are tilted to preserve the luminosity. The 
resulting primary vertex distribution, with cr^ < and > ax, is shown (the red 
ellipse). 

following formula: 



This gives, for all considered collider energies, a'^ ~ 3.6 /xm. This value is around 36 times 
greater than the spread axjV2 expected in case of collinear beams and comparable to the 
precision expected in the vertex position reconstruction. Therefore horizontal spread of 
the interaction point position cannot be neglected. 

For the analysis presented in this thesis the longitudinal size of the collision region 
is most important. As this is of the order of 100 /xm, we can expect that additional 
tracks and clusters due to overlaying events (resulting in additional vertexes, changed 
jet characteristics etc.) can influence the flavour-tagging algorithm and affect the event 
selection. Therefore, generation of all event samples used in the described analysis took 
into account the Gaussian smearing of primary vertex with aj, = 3.6 /xm, a'y = 11 nm and 
ct(, = 0.21 mm, and the beams crossing angle in horizontal plane, Oc = 34 mrad. 

3.4 The detector at TESLA 

The basic design of the TESLA detector for the Photon Collider is the same as for the 
TESLA mode. However, some modiflcations are needed due to the more complicated beam 
delivery system, including optical system which guides the laser beams to the conversion 
point. To protect detector components against the high-intensity low-angle radiation^ 
the tungsten mask is placed between the beam system and the detector. In case of the 

"^Background arises from synchrotron radiation and from upstream or downstream sources of 7 , e’*', 
e~ and n. 
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Figure 3.7: A scheme of one quadrant of the TESLA Detector ra. Dimensions are in 
mm. 


collider the opening angle of the mask is Omask = 83 mrad. Particles produced at 
smaller angles will not enter the detector. In case of the Photon Collider the value of 
Gmask = 130 mrad (7.5°) has been chosen as more space is required for optical system 
and beams removal. This results in moderate loss of hermeticity in comparison with the 
e’''e“-detector. Moreover, in case of the e’''e“-detector two forward calorimeters (Low 
Angle Tagger and Low Angle Calorimeter) are foreseen which together cover the region 
down to around 5 mrad. These components will not be installed in the Photon Collider 
option. 

In the following, main components of the detector for the Photon Collider are de¬ 
scribed. The description is based on the TESLA Technical Design Report (TDR) [201 
and the manual for the fast-simulation program Simdet EH. Eor many detector com¬ 
ponents different choices of technology and/or design were considered in the TDR. We 
discuss only these solutions which have been implemented in the SlMDET program and 
can be used to simulate the response of the detector. Because all proposed designs are 
expected to fulfill performance standards described in the TDR one can assume that our 
physical results should not worsen if alternative designs of the considered subdetectors 
are included in the final project. 
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The schematic view of the TESLA detector is shown in Fig. Id.71 The detector closest 
to the interaction point is a multi-layer microvertex detector with a total length of around 
30 cm. Currently at least two technologies are considered for this detector: charge-coupled 
devices (CCD) and active pixel sensors (APS). In this work the CCD-based option is used. 
It has well-defined geometry, material budget and the highest established performance in 
terms of precision over a wide range of incident angles (for devices of the dimensions 
needed for this application, i.e. tens of cm^), and only for this version the fast-simulation 
program provides a parametrized track covariance matrix which is crucial for the realistic 
flavour tagging simulation. With this design precision of the position measurement of 3.5 
yUm can be achieved. ^ To obtain a good reconstruction efficiency at least three detector 
layers are proposed so that, together with the silicon tracking subsystem (SIT), at least 
five silicon layers inside the TPC are available. 

In addition to the vertex detector a tracker system consists of intermediate silicon 
tracking detectors (SIT), a large Time Projection Chamber (TPC) and forward cham¬ 
bers. Silicon tracking subsystem includes cylinders in the barrel (SIT) and disks in the 
forward region (FTD). In the barrel region two layers of silicon strip detectors cover the 
region down to 6 = 25°. One of the cylinders, at r = 16 cm, improves the track recon¬ 
struction efficiency mostly for long-lived particles which decay outside the vertex detector. 
Three pixel and four strip silicon detectors with point resolutions of 10 /rm and 50 /rm, 
respectively are placed on each side, in the forward region (the endcaps). The main role 
of these detectors is to improve the momentum resolution for tracks by adding a few very 
precise space points at comparatively large distance from the primary interaction point, 
and to help the pattern recognition in linking the tracks found in the TPC with tracks 
found in the vertex detector. 

The central tracking system consists of two gas-filled chambers: a large volume time 
projection chamber (TPC) and a precise forward tracking chamber (FCH) located between 
the TPC endplate and the endcap calorimeter. The TPC, with 200 readout points in 
the radial direction (r = 32-170 cm), provides a very precise measurement of a track 
curvature, which is used in the determination of particle momentum. Because of the 
high magnetic field of 4 T the minimal transverse momenta of a particle required to 
enter and traverse the TPC are around 200 MeV and 1 GeV, respectively, if the particle 
charge is equal to the electron charge. Precise measurement of the specific energy loss 
in the TPC can be also used for particle identification. For tracks traversing the TPC 
at large polar angles the expected errors on the transverse momentum and the energy 
loss measurements are ail/px) ~ 0.02%/GeV and o'(^)/^ ~ 5%, respectively. For 
example, in case of a charged particle with energy of 20 GeV at the polar angle of 90° the 

®CCD vertex design implemented in Simdet assumes the radius of the innermost layer of 1.5 cm, 
which is the optimum choice for e“''e“ based on the background considerations. In case of the Photon 
Collider the inner radius of the vertex detector will probably have to be increased to about 2 cm. 
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momentum resolution is about 80 MeV. From ionisation losses the separation of kaons and 
pions should be possible in the momentum range from 2 to 20 GeV. Electron identification 
will be improved compared to what can be done with calorimeters alone, especially for 
low momenta (p < 3 GeV) where calorimetric identification is difficult. 

Overall tracking system performance, when the track parameters are determined from 
combining vertex detector, SIT and TPC measurements, shows a very high precision of 
transverse momentum determination a{l/pT) < 0.005%/GeV if systematic errors < 10 
/im for point position measurements are achieved. It is worth noticing that the overall 
momentum-resolution has been improved by about 30% by adding a cylindrical silicon 
detector (SIT) inside the TPG, i.e. at r = 30 cm. 

A tracking electromagnetic calorimeter (ECAL), build of tungsten absorber plates and 
thin silicon sensors, is placed behind the TPG. The expected energy resolution oeIE is 
around 11-14 %/a/e 7G^, depending slightly on the energy. The project assumes very 
high 3D granularity of this detector, allowing measurement of the particle momentum 
direction. A hadronic calorimeter (HGAL) is an iron/scintillating tile calorimeter with 
fine transverse and longitudinal segmentations. The energy resolution for single hadrons, 
estimated from simulation of hadronic showers in both calorimeters (HCAL+ECAL), is® 
geIE = 35%/v/^/GeV © 3%. 

A large superconducting coil, 6 m in diameter, produces a field of 4 T with very high 
uniformity (< 10“^). The coil is placed behind calorimeters to preserve high precision of 
energy measurement, reducing the amount of inactive material in front of the calorimeters. 
The iron return yoke serves also as a muon “separator”, absorbing other particles escaping 
from the HGAL. 

Eor muon chambers, which are placed inside the yoke and behind it, resistive plate 
chamber (RPG) technology is considered. Although the basic task for a muon detector is 
to identify muons, it is also possible to use the inner muon chambers as the “tail catcher”, 
i. e. to detect hadronic cascades which are not fully contained in the hadronic calorimeter. 
Eull efficiency of muon identification is reached for muons with energy above 5 GeV. 

The total length and the diameter of the detector will be around 15 m each. In gen¬ 
eral, the detector is designed to measure particles properties with a very high accuracy 
in the collision energy range from about 90 GeV up to 1 TeV. Electrons below 150 GeV, 
muons and charged hadrons are best measured in the tracking detectors. Electrons above 
150 GeV and photons by the electromagnetic calorimeter and neutral long-living hadrons 
by the combined response of the electromagnetic and hadronic calorimeters. In the event 
reconstruction the so-called energy-flow technique will be used which combines the infor¬ 
mation from tracking system and calorimeter to obtain the optimal estimate of the energy 
flow of produced particles and of the original parton four-momenta. Eor the energy-flow 

®The operator 0 means “adding in quadrature”: a(Bb = \/aE+lE 
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objects an average energy resolution of (Te/E ~ 30%/y/i?/GeV is expected. 

Due to the sparse beam structure (a long time interval of 199 ms between two bunch 
trains, a separation of two bunches inside a train by 337 ns, a train length of 950 /rs) no 
hardware trigger is foreseen. A total data volume of roughly 300 TByte per year will be 
stored for physical analyses. 

3.4.1 Simulation setup 

The fast simulation program for the TESLA detector, Simdet version 4.01 inn , was 
used to model the detector performance. All detector components are implemented in 
the program according to the TESLA TDK. Parametrizations based on the full simula¬ 
tion of detector performance with the Brahms program m are used to describe energy 
and angular resolutions. The track reconstruction efficiency and charge misinterpretation 
probability are momentum dependent. An energy-flow algorithm is used to link infor¬ 
mation from tracking system and calorimeters. In the first stage energy deposits in the 
calorimeters are joined into clusters. Then energy flow objects are defined by linking 
clusters with tracks reconstructed in the tracking system^. Eor the reconstructed particle 
tracks a track covariance matrix is calculated in the base: xy, z, 9, 0, 1/pt- 

Because two forward calorimeters, Low Angle Tagger and Low Angle Calorimeter, 
cannot be installed in the detector at the PC, they are not used in our simulation setup. 
Also the information about energy loss measured in the TPC, which is not properly simu¬ 
lated yet, is not used for particle identification. Instead, an appropriate misidentiflcation 
probability is assumed for each particle species. 

Within the current SlMDET version it is not yet possible to set a wider opening angle 
of the forward mask as required for the Photon Collider. To take the modified mask 
setup into account all generator-level particles are removed from the event record, before 
entering the detector simulation, if their polar angle is less than 9mask = 130 mrad. 


^In the current SlMDET version the idealised pattern recognition is still used, i.e. clusters are linked 
with tracks relying on the information about originally generated particles. 
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Chapter 4 

Signal and background 


In this Chapter the signal and background processes, and methods used in their simulation 
are described. The signal of the Higgs-boson production considered in this thesis is the 
process 77 —> higgs bb, whereas the main background processes are 77 —> bb and 
77 —cc. The bb pair production is an irreducible background which can be suppressed 
only by kinematic cuts. The cc pair production contributes to the background due to the 
finite probability of being tagged as the bb pair production. Background contributions 
from light quark and tan pair production, 77 — qq/T^T~ (q = u, d, s), are also considered. 
The process 77 —hadrons is described in detail due to its contribution to overlaying 
events. For heavy Higgs bosons, M^iggs > 160 GeV, also the process 77 —> W^W~ is 
taken into account. The underlying statistical principles used to describe the possibility 
of having more than one collision in single bunch crossing, and used in generation of the 
overlaying events are described in Appendix El 


4.1 Signal processes 

As the Higgs boson is a spinless particle, the distribution of its two-body decay products 
is uniform in the three-dimensional space (in the Higgs-boson rest frame). In spherical 
coordinates the distribution is uniform in cos 6 ** and 0 *, where 6* and 0 * are polar and 
azimuthal angle, respectively. In the Photon Collider higgs will be produced in collisions of 
photons which will have, in general, different energies. Thus, the center of mass system will 
be boosted with respect to the laboratory frame, resulting in a non-uniform distribution in 
cos 6 ', where 9 is the polar angle in the laboratory frame. However, as already mentioned 
in section mu energy of the electron beam is assumed to be tuned for the highest 
resonance production rate. It means that most collisions will involve two photons with 
high and similar energies as in all considered cases the resonance is narrow (the total width, 
^higgs, IS much Smaller than the width of the H0.^-distribution in the high energy part). 
Consequently, the Higgs boson will have small longitudinal momentum in comparison to 
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Symbol 

/i [GeV] 

M 2 [GeV] 

Aj [GeV] 

Mj [GeV] 

/ 

200 

200 

1500 

1000 

II 

-150 

200 

1500 

1000 

III 

-200 

200 

1500 

1000 

IV 

300 

200 

2450 

1000 


Table 4.1: MSSM parameter sets used in the described analysis. 

the mass, and its decays will be nearly isotropic also in the laboratory frame. 

Total widths and branching ratios of the Higgs bosons were calcnlated with the pro¬ 
gram Hdecay PI (version 3.0), where higher order QCD corrections are included. The 
mass of the top quark equal to 174 GeV was assumed. The contributions from the decay 
higgs gg* gbb were not added to the branching ratio BR{higgs bb) as the kine- 
matical characteristic of such events is different from the direct decay to the bb pair.^ For 
the considered mass range between 120 and 160 GeV the total width of the SM Higgs 
boson increases from abont 3.6 to 77 MeV, and its branching ratios BR(h — 77 ) and 
BR(h. ^ bb) decrease from 0.22% to 0.06% and from 68 % to 4%, respectively, in the mass 

Event generation for Higgs-boson prodnction process was done with the Pythia pro¬ 
gram (m. A parton shower algorithm implemented in Pythia was used to generate the 
final-state partons. The fragmentation into hadrons was also performed using the Pythia 
program, both for Higgs-boson prodnction and for all backgronnd event samples. 

In case of the MSSM Higgs-boson production the analysis has been developed assuming 
MSSM parameters similar to these nsed in [2lj, i.e. tan/3 = 7, /i = ±200 GeV and 
M 2 = 200 GeV, taking into acconnt decays to and loops of snpersymmetric particles. The 
parameter value tan /3 = 7 is used in the event generation and obtained results are rescaled 
to the parameter range tan/3 = 3 — 20. In the following, parameter sets from |2l| will be 
denoted as / and III, see Tab. mu Only the valnes of trilinear conplings are changed 
(from Aj = 0 to Aj = 1500 GeV), so that the mass of the lightest Higgs boson, instead 
of being around 105 GeV (for tan/3 = 4 and Ma = 300 GeV) is above the current lower 
limit for the SM Higgs boson, > 114.4 GeV. Resnlts for heavy nentral Higgs bosons 
are the same as for parameter sets proposed in 121 . The intermediate scenario II with 
fi = —150 GeV is also considered. For comparison with predictions presented by LHC 
experiments, the scenario IV nsed in is also inclnded. In all cases the common sfermion 
mass eqnal to 1 TeV was assumed. We have checked that all parameter sets imply masses 
of neutralinos, charginos, sleptons and squarks higher than current experimantal limits. 

^Inclusion of the decay higgs gg* gbb does not change the results of this analysis as BR(/i ^ 
bb) increases only by around 1% for the SM case, and for considered MSSM parameters BR{A/H —s- 
gg)/BR{A/H bb) < 1%. However, if events with only one 6-tagged jet were accepted, one would have 
to use inclusive branching ratio BR(h —> bb + ...), measured at the e+e” LC, to obtain final results for 
r{h 77). 
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Figure 4.1: The total cross section for process 77 —> h —> 66 for set I of the MSSM 
parameters values (see text) with Ma = 300 GeV and tan/? = 7. The signal only, as, and 
interference part, < 7 /, are shown separately. As the other neutral Higgs bosons are narrow 
and heavy, they do not influence the h cross section. 

For the wide range of parameter values the heavy neutral Higgs bosons, A and H, are 
nearly mass degenerate. The mass difference Mh — Ma decreases with increasing tan/? 
and Ma and is similar for all considered parameter sets. For Ma = 200 GeV the mass 
difference decreases from Mh — Ma ~ 12 GeV for tan/? = 3 to 0.7 GeV for tan/? = 15. 
If Ma = 350 GeV, the corresponding values are 6 GeV and 0.3 GeV. The mass difference 
is larger or comparable to the total widths of A and H which vary between 50 MeV and 
4 GeV. The branching ratios relevant for this study change between 3% and 90% for 
BR{A/H —> 66 ), and from 2 ■ 10“^ to 9 ■ 10“® for BR(A/H —> 77 ). 

Processes 77 ^ 77 —66 and 77 —A —66 do not interfere^ {e.g. see [22] eq. 3.15), 
and h — H interference is negligible due to the large difference in masses and relatively 
small widths of the bosons. Therefore total 77 —> A, 77 —66 production rate is equal to 
the sum of both contributions. 

However, for the complete description of the Higgs boson production we also have to 

^Although in general case two processes with A and ff in the intermediate state (different CP quantum 
numbers) can interfere. For example, there is interference between chargino production processes in 77 
collisions via A and H: ^ A, H ^ xtxj' R-9- see [22] eq. 3.59). 
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Figure 4.2: The total cross section for process 'f'y ^ A, H ^ bb for set I of the MSSM 
parameters values (see text) with Ma_ = 300 GeV and tan/? = 7. The signal only, as, 
and interference part, ai, are shown separately. As the light Higgs boson, h, has mass of 
aronnd 120 GeV and is narrow, it does not inflnence the A, H cross section. 


consider the interference between 77 —higgs bb and non-resonant 77 —bb production 
processes. The LO interference terms for h and H/A production are shown in Fig. 14.11 
and 14.21 together with the signal cross sections 77 —> h —> 66 and 77 —> A, if ^ 66 . For 
all considered cases these terms are proportional to the real part of the propagator: 

^ m2 M2(s-M2) 

s - M2 + 5MF (s - M2)2 + (MF)2 

As on average (the integral over = ^/s) this expression is near to zero for small F 
(of the order of F/M), the interference contribution can be safely neglected . However, 
as discussed in the next section, the NLO corrections substantially modify predictions for 
heavy quark production. The NLO corrections for interference term were calculated in 
[ 22 ] . As described in j 22 j, after selection cnts the interference term was below the level 
of 10“^ of the signal. Although selection cuts were different than in our analysis, we 
can infer the total correction factor. The cuts decrease interference contribution by the 
order of magnitnde and the signal rate by about 50%. Thus, one can estimate that the 


34 












interference part contributes no more than 1% of the total signal cross section. Because 
this is smaller than other uncertainties, we neglect the interference term in this analysis. 

4.2 Heavy quark production background 

The main background for the considered signal process, 77 —> higgs —> bb, is the heavy 
quark-pair production. An irreducible background consists of events with the bb{g) final 
state, resulting from ’direct’ nonresonant bb production, 77 — bb. In LO approximation 
the cross section for 7,^ = 0 is suppressed and the dominant contribution is due to the 
\Jz\ = 2 state. This is very fortunate as the 77 -luminosity spectrum is optimized to give 
highest = 0 luminosity and the \Jz\ = 2 component is small in the fii^'^'s-production 
region. Unfortunately, NLO corrections compensate partially the mg/s-suppression and, 
after taking into account luminosity spectra, both contributions (for Jz = 0 and | Jz| = 2 ) 
become comparable. The extensive comparison of NLO and LO results can be found for 
example in our work 123 

The other processes 77 —> qq{g), where q = u,d,s,c, contribute to the reducible 
background. However, one has to consider these processes due to the non-zero probability 
of wrong fiavour assignment in reconstruction (impurity of flavour-tagging). Events with 
cc{g) in the final state have the highest mistagging probability. In comparison to the 
77 —> bb process there is an enhancement factor of (cc/ef,)"^ = 16 in the 77 —*• cc cross 
section. It turns out that after fiavour tagging both processes give similar contribution to 
the background. 

The background events due to processes 77 —> bb{g), cc{g) were generated using the 
program written by G. Jikia [23], where a complete NLO QCD calculation for the produc¬ 
tion of massive quarks is performed in the massive-quark scheme. The program includes 
exact one-loop QCD corrections to the lowest order processes 77 —*• bb, cc na. and the 
non-Sudakov form factor in the double-logarithmic approximation, calculated up to four 
loops | 2 I|. Events generated with NLO QCD program were transfered to Pythia program 
for hadronisation. To avoid double-counting of corrections due to real gluon emission, the 
parton shower algorithm was not applied. However, to estimate the influence of higher 
order corrections on the event selection efficiency we also prepared dedicated samples of 
77 —> bb{g), cc{g) events with parton shower included. Results of this comparison are 
presented in Appendix [B] 

4.3 Other background processes 

In cases of the SM Higgs-boson production for Mh = 150 and 160 GeV, and in the 
analysis of heavy neutral Higgs-bosons in the MSSM also the pair production of lU bosons. 
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77 IV~^IV~, is considered as a possible background. The cross section for this process is 
very high for large ITiy-y, and it can contribute to the background if the event is clustered 
to two or three jets and at least one of these jets is 6 -tagged. This can be the case if 
two jets from hadronic IV decays are merged together by the jet-clustering algorithm, 
or if some jets are ignored in the analysis because they are too close to the beam pipe. 
For generation of 77 —> W~^W~ events the Pythia program is used. However, as only 
unpolarized cross section for this process is implemented in Pythia, we use polarized 
differential cross section formulae from |S 3 | to obtain correct distributions for Jz = 0 and 
\Jz\ =2 contributions. 


As there is non-zero probability of mistagging a light-quark jet as a &-jet, the process 
77 ^ QQi where q = u,d, s, is also taken into account as a possible background. Due to 
the strong dependence of the cross section for this process on the fermion charge (a ~ e^), 
the 77 —> uu contribution dominates. The event generation is performed with PYTHIA 
using unpolarized LO cross section. We known that for Jz = 0 the LO cross section for the 
process 77 —qq is equal to zero for massless quarks. By convoluting the cross section for 
IJ 2 I = 2 with the total luminosity spectrum (modulo factor 2 ) we overestimate the light 
quark production background. Comparing results for 77 —^ cc{g) and 77 —> uu we have 
determined that the number of events with light quark-pair production is overestimated by 
a factor of about 2.6 for our SM analysis, and by a factor of about 4 for our MSSM analysis. 
We do not apply any corrections to decrease this effect. In the analysis of SM Higgs-boson 
production the light-quark contribution is negligible and does not change the results. In 
case of MSSM Higgs-boson production we use the overestimated contribution of 77 ^ qq, 
q = u,d, s, events to effectively take into account the contribution of 77 —hadrons events 
(without directx direct interactions) from which no generated events passed all selection 
cuts. We estimated that in the mass window optimal for the cross-section measurement 
(see Chapters El and Ej) the contribution of 77 —hadrons events would correspond to 
about 50% of the (overestimated) light-quark contribution. 


A large number of r-pair production events will also be observed in the Photon Collider 
(the cross section about 1.7 times higher than the LO cross section for 77 —>■ cc). Thus, 
even small r mistagging probability could in principle significantly influence our results. 
Fortunately, the Jz = 0 contribution is negligible in this case due to the m^/s-suppression 
(QED higher order corrections are small). In addition, considerable amount of energy 
is carried out by neutrinos. Consequently, the number of r+r“ events reconstructed as 
two-jet events with high invariant mass should be very small. For this study 77 —r+r“ 
production events were generated with Pythia using polarized cross section. 
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4.4 Overlaying events 77 ^ hadrons 


Because of the large cross section and huge 77-luminosity at low from one to two 
77 —> hadrons events'^ are expected at the TESLA Photon Collider per bunch crossing. 
These events hardly contribnte to the background on their own. However, they can have 
a great impact on the reconstruction of other events produced in the same bunch crossing, 
by changing their kinematical and topological characteristics. 

We generate 77 ^ hadrons events with Pythia, using the Inminosity spectra from 
a full simulation of the photon-photon collisions IHEl, rescaled to the chosen beam en¬ 
ergy. For each considered e“e“ energy, an average number of the 77 —> hadrons 

events per bnnch crossing is calculated. Then, for every signal or backgronnd event, the 
77 — hadrons events are overlaid (added to the event record) according to the Poisson 
distribution. In Appendix El principles of event generation and approximations used are 
discussed in detail. 

The 77 — hadrons processes are classified according to the type of photon(s) inter¬ 
action. If both photons interact as point-like, as described by QED, then we call the 
process directx direct. However, if one or two photons interact as a hadronic state (vector 
meson or qnantnm finctnation with gluons), then they are denoted as hadron-likex direct 
or hadron-likex hadron-like, respectively. As shown in Fig. 14.31 the biggest contribntion to 
the cross section is due to processes with hadron-likex hadron-like photons. Fortunately, 
the 77 — hadrons cross section is very forward-peaked as seen in Fig. 14.41 A cnt on 
the polar angle of tracks and clnsters measured in the detector should greatly reduce 
contribution of particles from 77 —> hadrons processes to selected events. Events with 
77 —> hadrons interactions only, withont directx direct interactions, are simnlated as well. 
Snch events can mimic signal if two or more of them are overlaid. As the generation of 
high-energy 77 — hadrons events with significant transverse energy is very inefficient, 
these events are included only in the SM analysis. In case of MSSM analysis we have 
estimated the contribntion of 77 ^ hadrons events, and inclnde it effectively as the part 
of light quark-pair production (see previous Section). 

For more details concerning 77 —> hadrons overlaying events and their influence on 
the reconstrnction see Appendix O The package Orlop, created for inclnding the 
77 —> hadrons events by the generation of hard 77 scattering processes, is described 
in Appendix ini 


^For technical reasons we consider only photon-photon events with > 4 GeV. However, events 
with lower fFy.y are mostly produced with high boost and particles going at very small angles do not 
enter the detector. For further detailed discussion see Appendix O 
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PYTHIA 6.214: Total cross section yy ^ hadrons 



Figure 4.3: Cross sections for the process 77 —hadrons obtained with the Pythia 
program. Contributions from various event classes are indicated. Also the magnitude of 
the cross section for elastic and diffractive processes is indicated. 
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Transverse energy per bunch crossing 
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Figure 4.4: Angular distributions of the transverse energy flow, E^, for overlaying 

events 77 —hadrons per bunch crossing. Various components and their total contribu¬ 
tions (numbers in parentheses) are indicated. Pythia generation results with luminosity 
spectrum for = 210.5 GeV. 
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Chapter 5 

Standard Model Higgs-boson 
production 


In this Chapter the measurement of the Standard Model Higgs-boson production cross 
section at the TESLA Photon Collider is discussed. Following steps of the analysis are 
described: selection of energy-flow objects, jet reconstruction, kinematical and topolog¬ 
ical selection cuts optimized for cross section measurement, and the role of 6 -tagging in 
selection of signal events. To simplify the description, the analysis is presented in detail 
for the Higgs-boson mass Mh = 120 GeV. For other considered masses of the Higgs boson, 
Mh = 130, 140, 150 and 160 GeV, the same procedure was performed with indepen¬ 
dent optimization of selection thresholds. The cuts dedicated to suppress 77 —> W^W~ 
background, which is not relevant for lower Higgs-boson masses, are described for Mh = 
160 GeV. Expected precisions of the measurement obtained in this analysis are compared 
with results of our earlier works in which some of experimental aspects and background 
contributions considered here were not yet taken into account. 


5.1 Preselection of energy-flow objects and jet recon¬ 
struction 

In the energy range = 210-260 GeV about one 77 —> hadrons event takes place on 
average at each bunch crossing. The contribution from these overlaying events is expected 
to affect observed particle and energy flow mainly at low polar angles (see Section HI. 
Therefore, we introduce an angle 9tc defining the region strongly contaminated by this 
contribution; tracks and clusters with polar angle less than 9tc are not taken into account 
when applying energy-flow algorithm. In spite of that energy-flow objects with polar angle 
less than 9tc can still be formed; they are also ignored in further steps of analysis. A few 
values of 9tc were considered in the analysis as discussed in detail in Appendix O We 
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decided to use the value 6tc = 0.85 as with this choice almost the whole contribution 
from hadron-like photon interactions is suppressed and distributions of jet transverse 
momentum and jet mass are similar to those obtained withont overlaying events and 
without 9tc cut. It was checked that 9tc = 0.85 results also in the best final cross 
section measurement precision. 

For the signal process considered in this analysis we expect that the produced partonic 
state is well reproduced by jets reconstructed from energy-fiow objects. In the presented 
study jets are reconstructed using the Durham algorithm lEH where the distance measure 
between two jets, i and j, is defined as 

yij = - cos9ij)/E‘^; 

Ei and Ej are energies of jets, % is the relative angle between jets and E is the total 
energy measured in the detector. The list of energy-fiow objects reconstructed in the 
detector is used as the input to the algorithm, assuming that each energy-fiow object is 
a jet. In following steps a pair of jets which has the smallest value of yij is searched for 
and these two jets are merged into one jet. The algorithm terminates when all possible 
values of y^j are greater than the value of the cut-off parameter, ycut- The choice of the 
distance measure yij and of the parameter ycut value used in this analysis is based on 
the approach adopted in the NLO QCD 121 program which is used for generation of 
backgronnd events 77 ^ QQ{g)- In this program the real gluon emission is considered 
only for ygg > 0.01. Soft gluon emissions, i.e. emissions with ygg < 0.01, are absorbed in 
the cross-section calculation for QQ final state. For consistency with this approach jets 
have to be reconstructed with ycut > 0 . 01 , as for lower ycut values additional jets expected 
from soft-gluons emission would not be described by the generator. Moreover, the distance 
measure used in the NLO generator is calculated using true values of kinematic variables 
and is inversely proportional to the 77 invariant mass squared, oc whereas the 

visible energy is used in the jet reconstruction, yij oc l/E"^. In the significant fraction 
of events we expect that due to detector acceptance E'^ < s^-^. Therefore, the value 
ycut = 0.02, two times larger than the one used in generator, has been chosen. With this 
value, reconstructed jets can be relatively wide. For example, two perpendicular jets will 
be joined together if one of them has energy of 12 GeV or below (assuming E 120 GeV, 
most probable value for Mh = 120 GeV). 

As the NLO QCD generator does not include additional gluon emissions due to higher 
order corrections, we study the resulting systematic uncertainty of the result, by applying 
the parton shower algorithm to the NLO heavy quark background events. Although 
some gluon contributions are double counted in such procedure, it allows us to determine 
the sensitivity of the analysis to the higher order corrections. Results are presented in 
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Appendix |Bl 

To correct for the non-zero beam crossing angle, all reconstructed jets are transformed 
from the laboratory frame to the frame moving with the speed factor /? = sin ^ ^ = 

0.017 in the x direction, where Oc is the beam crossing angle. After this correction the 
average value of the measured transverse momentum in the horizontal direction, P^, is 
zero. 


5.2 Kinematical and topological cuts 

The first cut applied after detector simulation is introduced to exclude possible influence 
of the cut IT^™, the lower limit on the 77 invariant mass, used in the event generation. 
Therefore, the condition Wrec > 1.21T™“ is imposed for all considered events, where IVrec 
is the total reconstructed invariant mass of the event (calculated from all energy-flow 
objects above 0tc)- 

Higgs-boson decay events are expected to consist mainly of two 6 -tagged jets with 
large transverse momentum and nearly isotropic distribution of the jet directions. The 
significant number of events (~ 25%) contains the third jet due to the real gluon emissions 
which are approximated in this analysis by the parton shower algorithm, as implemented 
in the Pythia. 

The following cuts are used to select properly reconstructed bb events coming from 
Higgs decay. 

1. Number of selected jets should be 2 or 3. In addition to two 6 -quark jets we allow 
for one additional jet from hard gluon emission. The signal-to-background ratio 
is similar for both jet multiplicities. Moreover, the NLO QCD generator used for 
heavy-quark background generation does not include resummation of the so-called 
Sudakov logarithms which would be relevant if 2- and 3-jet events classes were 
considered separately. 

2. The condition | cos< Cq is imposed for all jets in the event where is the jet 
polar angle, i.e. the angle between the jet axis and the beam line. This cut should 
improve signal-to-background ratio as the signal is almost uniform in cos^, while 
the background is peaked at | cos^l = 1 . 

3. Since the Higgs bosons are expected to be produced almost at rest, the ratio of the 
total longitudinal momentum calculated from all jets in the event, P^, to the total 
energy, P, should fulfill condition |P^|/P < Cp^. 

To determine the cut parameter values Cq and Cp^ the corresponding distributions of 
the signal and heavy quark background events were compared (other background contri¬ 
butions were not considered at this stage). After cut 1 the optimal value of parameter Cq 
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Figure 5.1: Distributions of | cos|(maximal value of Icos^jetl over all jets) for 
signal and background events. For background only 77 —^ bb{g) events are shown. The 
signal measurement precision /\a/a of about 7% is obtained for | cos 6 'jet|™'^'" < Cq = 
0.725. 


(cnt 2 ) is fonnd as the one which minimizes the estimated statistical nncertainty of the 
measurement: 

Aq-(77 -^h^bb) _ y/fis + Fb 
< 7(77 ^ h ^ bb) jjLs 

where ps and fis are the nnmbers of expected signal and backgronnd events after the cnt, 
respectively. With optimized cnt 2 the same procedure is repeated for parameter Cp^ (cut 
3). The expected event distributions for | cos9jet\'^^^ (the maximum value of | cos 9jet\ over 
all jets in the event) and \Pz\/E are shown in Fig. 15.1 l and 15.21 respectively. For simplicity 
only 77 —> bb{g) backgronnd contribntion is shown. The 77 —> cc{g) contribntion, which 
is around 16 times larger, has a very similar shape. Both background contributions are 
taken into account in cut optimization. For Mh = 120 GeV the optimized cut values are 
Cq = 0.725 and Cp^ = 0.1, as indicated in the fignres (vertical arrows). The measnrement 
precision is estimated to be around 7% and 5% after the | cos^F*] cut and after the \Pz\/E 
cut, respectively. Angular cuts used in the event selection procedure are compared in Fig. 
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Figure 5.2: Distributions of \Pz\/E for signal and background events. For background 
only 77 ^ bb{g) events are shown. The signal measurement precision lS.a ja of about 5% 
is obtained for \Pz\/E < Cp^ =0.1. 


5.3 6-tagging algorithm 

For 6-tagging the Zvtop-B-Hadron-Tagger package prepared for the TESLA project 
was used |HK1 inni IH7j . The flavour tagging algorithm is based primarily on ZVTOP, the 
topological vertex finding procedure developed at SLD IHHl- In addition to the ZVTOP 
results, a one-prong charm tag 1^11 and an impact parameter joint probability tag |fi9| 
outputs are used to train a neural net. Following parameters are given as an input to the 
neural-network algorithm (for all tracks or all vertices): 

1. Impact parameters in r —0 and r — z. Impact parameter in the r —0 plane is defined 
as the minimal distance between the track trajectory and the beam axis; impact 
parameter in r—z plane is defined as the distance between the reconstructed primary 
vertex position and the point on the beam axis nearest to the track trajectory. 

2. Significance of the track impact parameters - the ratio of the impact parameters to 
their estimated errors. 

3. Vertex decay length - the distance between the primary vertex and the secondary 
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Figure 5.3: Comparison of the angular cuts used in the selection procedure of the 

Higgs-boson production events for = 120 GeV. 

or tertiary vertex. 

4. Vertex decay length significance - the ratio of the vertex decay length to its mea¬ 
surement error. 

5. Pi-corrected mass of the secondary vertex - the invariant mass of particles coming 
from the vertex. As only charged particles (tracks) are considered, the correction 
for neutral particles is applied. The correction is based on the assumption that the 
total momentum of all particles coming from the secondary vertex must be parallel 
to the vector between the primary and secondary vertex positions. 

6 . Vertex momentum - the total momentum of all tracks belonging to the vertex. 

7. Secondary vertex track multiplicity. 

8 . Secondary vertex probability - the probability that all tracks assigned by the ZV- 
TOP algorithm to the secondary vertex belong to this one vertex. 

The neural-network algorithm was trained on the Z decays. For each jet the routine 
returns a “ 6 -tag” value - the number between 0 and 1 corresponding to “ 6 -jet” likelihood. 

In order to optimize the signal cross-section measurement, the two-dimensional cut on 
6 -tag values is used. In the signal events two jets with the highest transverse momentum 
are most likely to originate from 6 quarks. Therefore, all jets in the event are sorted 
according to the value of their transverse momentum. The distribution of 6 -tag values for 
2 and 3-jet events is considered in the plane 6 -tag(jeti)( 8 ) 6 -tag(jet 2 ) where indices 1 and 
2 correspond to two jets with the highest transverse momenta. The two-dimensional 
distributions of 6 -tag values for the signal, 77 ^ 6 , —> 66 , and for the background, 
77 —> 66 (p), events are shown in Fig. 15.41 The corresponding distributions for other 
considered background contributions, 77 ^ cc(p) and 77 —qq {q = u,d, s), are shown 
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Figure 5.4: Distributions of 77 —> /z. ^ 66 (left) and 77 ^ bb{g) (right) events in the 
plane btag(jeti) (81 btag(jet 2 ). 


in Fig. 15.51 Events considered in the 6 -tagging studies fulfill fore-mentioned, optimized 
selection cuts and an additional cut Wrec >0.7 Mh which removes low-mass events not 
relevant for the final result (this cut is used only for tagging optimization). As expected, 
for processes 77 —> 6 . ^ 66 and 77 ^ bb{g) most events populate the regions with high 
6 -tag values (Fig. 15.411 . whereas most 77 —> cc{g) and 77 — qq events have small 6 - 
tag values iFig. I5.5ll . Nevertheless, significant fraction of 77 —> cc{g) events populates 
the region of high 6 -tag values, and the event distribution is more fiat than the one for 
77 QQ events. The optimal higgs-tagging cut is found by considering the value of 
the signal to background ratio S/B, where S and B denote the expected numbers of 
events for the signal 77 ^ 6 . —> 66 and for the sum of background contributions from 
processes 77 — QQ{g) {Q = c,b) and 77 — qq {q = u,d,s), respectively. Obtained S/B 
distribution in the 6 -tag(jeti)® 6 -tag(jet 2 ) plane for Higgs-boson production with Mh = 
120 GeV is shown in Fig. 15.61 The selection criteria which results in the best precision of 
the F(/z —> 77 )BR( 6 , — 66 ) measurement corresponds to S/B > 0.19 as indicated in the 
figure (stars). 

The obtained efficiencies for tagging higgs events, 66 background events, and the prob¬ 
abilities for mistagging of the cc and qq {q = u, d, s) events are Eh = 58%, Ebb = 50%, 
Ecc = 2.2% and Ends = 0.16%, respectively. Similar efficiencies were obtained for other 
considered electron-beam energies. We call this procedure ’hz^'^'s-tagging’ as the efficiency 
for tagging signal events. Eh, is significantly higher than the efficiency for tagging bb{g) 
background events. Ebb- This is because large fraction of background events is recon¬ 
structed as 3-jet events (LO contribution is suppressed for = 0) in which the gluon 
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Figure 5.5: Distributions of 77 —*• cc{g) (left) and 77 —> qq, q = u,d,s, (right) events in 
the plane btag(jeti) ® btag(jet 2 ). 

jet is often one of the two jets with highest transverse momenta. In the earlier analyses 
123 EH a fixed 66 -tagging efficiency, Sh = Ebb = 70%, and a fixed cc-mistagging efficiency, 
Ecc = 3.5%, were assnmed. Althongh the efficiencies resnlting from the optimized selection 
are much lower, signal to background ratio Eh/Ecc improves significantly. 

Particles from 77 — hadrons overlaying events can significantly change properties 
of the jet to which they are assigned by the jet clustering algorithm. For example, the 
invariant mass of the jet increases on average by 3 GeV, if the angnlar cut is not applied 
{i.e. parameter 9tc = 9mask', see Fig. IC.1 fl in Appendix EJ. Althongh the average in¬ 
variant mass of the jet after the cut corresponding to cos 9 tc = 0.85 is similar to the jet 
mass without overlaying-events contribution, the jet structure can still be affected by the 
remaining particles from 77 ^ hadrons interactions, and by rejection of some particles 
coming from the signal process. These effects influence also the flavour tagging algorithm 
and cause a significant change in the results of the 66 -tagging optimization. To quantify 
the infinence of overlaying events we repeated the optimization procednre, for prodnction 
of the SM Higgs boson with Mh = 120 GeV, withont overlaying events and with cos9tc = 
0.99. The resulting efficiencies, corresponding to optimal cut S/B > 0.16, are Eh = 71%, 
Ebb = 64%, Ecc = 2.9% and Euds = 0.11%. The corresponding selection region in the 
6 -tag(jeti)® 6 -tag(jet 2 ) plane is significantly wider than for the nominal analysis, but the 
cc background suppression factor Ebb/Ecc is similar. 

Although the selection region is smaller when overlaying events are taken into account, 
the efficiency for 77 — qq, q = u,d,s, is greater by about 50%. Jets coming from light- 
qnark prodnction can be significantly modified by 77 — hadrons events, and are more 
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Figure 5.6: The expected ratio of signal (77 ^ h ^ bb) to background (77 —> QQ{g), 
Q = c,b, and 77 —> qq, q = u,d, s) event distributions in the plane btag(jeti) ® btag(jet 2 ). 
The region which results in the best precision measurement for the cross-section measure¬ 
ment is indicated by stars. 


likely to be recognized as 6-jets. One of important reasons is that the primary vertex 
of the overlaying event is usually shifted with respect to the primary vertex of the hard 
interaction. If particles from both interactions are combined in the reconstructed hadronic 
jet, the vertex finding algorithm is likely to reconstruct two vertices treating one of them 
as a primary vertex, and the second one as the vertex resulting from b decay. The effect of 
tagging deterioration is even stronger for higher beam energies as more overlaying events 
per bunch crossing are produced due to higher luminosity and cross section. The influence 
of overlaying events on the event reconstruction is clearly seen also in case of the heavy 
MSSM Higgs-bosons production as will be discussed in Chapter El 

The influence of the optimized higgs-tagging criteria on the reconstructed invariant- 
mass, Wrec, distribution is shown in Fig. 15.71 for signal events ■yy ^ h ^ bb after all 
described selection cuts. These results were obtained without overlaying events 77 —> 
hadrons. The tail towards low masses is due to events with energetic neutrinos coming 
from semileptonic decays of D and B mesons (see for more details). Contribution 
of these events can be suppressed by an additional cut Pt/Et < 0.04, where Pt and Et 
are the absolute values of the total transverse momentum of an event, P^, and the total 
transverse energy, respectively. We see that the efficiency of 6-tagging is similar for events 
with and without energetic neutrinos as the algorithm does not significantly influence the 
shape of the distributions. 
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Figure 5.7: Distributions of the reconstructed invariant mass, Wrec selected 77 ^ 
h ^ bb events, for Mh = 120 GeV. Distributions obtained before and after applying the 
6 -taggingalgorithm, without and with the additional Pt/Et < 0.04 cut are compared. 
Overlaying events are not inclnded. 

5.4 Results 

The invariant-mass distribntions for signal events passing all optimized selection cnts, 
before and after taking into account the overlaying events 77 —hadrons are compared in 
Fig. 15.81 Tefti. Mass resolution, derived from the Gaussian fit in the region from /i— l.Scx 
to /r+1.3(7, is abont 4 and 6 GeV, respectively. The overlaying events and cnts snppressing 
their contribution significantly influence the mass reconstruction and result in an increase 
of distribution width by about 2 GeV, and in a shift of the mean value, p., by about 3 GeV. 
A drop in a selection efficiency, resnlting in the rednced nnmber of events expected after 
selection cuts (from about 6450 to 5530 events), is also observed. This is because tighter 
6 -tagging cuts have to be imposed to reduce influence of overlaying events. More events 
are also rejected by the cnt on the longitndinal momentnm. Some drop in the selection 
efficiency is also dne to the fact that the energy deposits from the 77 —hadrons processes, 
remaining after the 9tc cut, “shift” jets nearer to the beam axis and the event can be 
rejected by the jet-angle cnt. 

Distribntions of the reconstrncted invariant mass, Wrec, expected after applying all 
selection cuts and 6-tagging algorithm, for the signal (77 h ^ bb) and all considered 
background contributions are shown in Fig. 15.91 Beside the heavy quark production, 
77 —> bb{g) and 77 —> cc{g), all background contributions are small, of the order of 30%, 
15% and 5% of the signal for 77 — hadrons, 77 —> qq and 77 —>• r+r“, respectively. 

We assume that the number of observed Higgs-boson production events will be ex- 
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Figure 5.8: Reconstructed invariant-mass , Wrec, (left) and corrected invariant-mass 
, Wcorr, (right) distributions for selected 77 ^ h ^ 66 events, for Mh = 120 GeV. 
Distributions obtained without and with overlaying events (OE) are compared. Results 
for the mean and dispersion a from the Gaussian fit in the region from fi — l.Scx to 
^ + 1.3(j, are also shown. 


tracted by counting the number of 66 events in the mass window around the Higgs-boson 
mass peak, Nohs, and subtracting the expected contribution of background events, hb- 
The relative statistical error expected in the measurement of the Higgs-boson production 
cross section ( 7(77 —> 6 , —> 66 ), or of the partial width multiplied by the branching ratio 
r(6, —77)BR(6, — 66 ), can be estimated from the following formula: 


A(t( 77 —6, —66) A [r(6, ^ 77)BR(6, —> 66)] 




< 7(77 h ^ bb) 


r{h 77)BR(fi ^ 66) 




( 5 . 1 ) 


where /i^ is expected number of signal events, and the expected number of observed events 
is (Nobs) = + fJ'B- The accuracy obtained from the reconstructed invariant-mass dis¬ 

tribution for Mh = 120 GeV is equal to 2.3%. The mass window used to calculate the 
signal measurement precision is again optimized to obtain the lowest relative error. For 
Mh = 120 GeV the selected mass window is 100 to 127.5 GeV as indicated in Fig. 15.91 
(vertical arrows). The obtained result is consistent with results of our previous analyses 
mm which however did not take into account many aspects of the measurement con¬ 
sidered here. In the current analysis additional background contributions deteriorate the 
precision of the cross section measurement. However, the effect is partially compensated 
by the performed optimization of selection procedure. 

Significant part of the energy in the signal events can be lost due to escaping neutrinos. 
As shown in Fig. EH this effect worsens the mass resolution and, as a result, the cross 
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Figure 5.9: Distributions of the reconstructed invariant mass, Wrec, for selected bb events. 
Contributions of the signal, for Mh = 120 GeV, and of the background processes, i.e. 
77 QQ for (l = u,d, s, 77 — > , and 77 ^ hadrons (as a separate contribution with 

hadron-likexhadron-like interactions only, indicated as ’resolved’), are shown separately. 
Arrows indicate the mass window, 100 to 127.5 GeV, optimized for the measnrement of 
the T{h 77)BR(/i ^ bb), which leads to the statistical precision of 2.3%. 


section measurement precision. We have looked for a correction method which would 
improve the mass resolution without reducing the event statistics. Unfortunately, due to 
a large spread of the photon beam energy, no nnambignons constraint can be imposed 
on the reconstrncted longitndinal momentnm. We have considered fonr methods of 
correcting the measured invariant mass: 

1. The valne of the measnred transverse momentnm is added to the total energy and 
the transverse momentnm is balanced. This is eqnivalent to the assnmption that 
the missing transverse momentum is due to a single neutrino emitted perpendic¬ 
ularly to the beam line. The procedure is illustrated in Fig. 15.101 where the 2-jet 
77 —> /r —66 event is shown with energy and momentnm carried ont by two nen- 
trinos. As the total longitudinal momentum, Pl, is not constrained, the missing 
longitudinal momentum, cannot be estimated. However, the total transverse 

momentnm shonld be balanced. Thns, the missing transverse momentnm, 
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Figure 5.10: The 2-jet 77 —> h ^ 66 event in which part of energy is carried out by 
neutrinos. As the total longitudinal momentum, Pl, is unknown, the missing longitudi¬ 
nal momentum, cannot be estimated. However, the total transverse momentum 

should be balanced and the missing transverse momentum, can be attributed to 

the ’effective’ neutrino, z/g//. 


can be attributed to the ’effective’ neutrino. This correction was introduced in our 
earlier analysis 123 . 

2. The transverse momentum, px, of the jet with lowest px is increased by a value 
of the total missing transverse momentum and the longitudinal jet momentum is 
rescaled to preserve its original direction. This method is applied for 2-jet events 
and assumes that the missing Px is due to the single neutrino emitted along the jet 
with lower px- In general the total transverse momentum is still unbalanced after 
this correction. 

3. The transverse momentum of 2-jet event is balanced under assumption that the 
missing Px is due to the neutrino emitted under the polar angle equal to the polar 
angle of the jet with lowest px- 

4. The transverse momentum of 3-jet events is balanced by rescaling momenta of two 
jets. All combinations are checked and the most “reasonable” one is chosen, i.e. the 
one satisfying the requirement that each of two rescaling factors is greater than 1 
and less than 1.3. 

Approaches combining these methods were also taken into account. Moreover, we also 
considered the algorithms where the correction were limited to events with transverse 
momentum greater than some threshold value, which was varied to obtain the best mass 
resolution. Surprisingly, the first (and simplest) procedure proved to be the best one when 
applied to all events. Other correction methods introduce systematic bias in the corrected 
invariant mass distributions and sometimes result even in the deterioration of the final 
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Figure 5.11: As in Fig. 15.91 for the corrected invariant mass, Wcorr, distributions. 

Arrows indicate the mass window, 107.5 to 132.5 GeV, optimized for the measurement of 
the F(h ^ 77 )BR(h ^ bb), which leads to the statistical precision of 2.1%. 


cross-section measurement precision. Thus, the corrected reconstructed invariant mass 
used for the final analysis was defined as 123: 

+ 2Pt(E + Pt). (6.2) 

In Fig. 15.81 (right) the distributions of Wcorr for the selected signal events, without 
and with overlaying events, are presented. The tail of events with invariant masses below 
~ 110 GeV is mnch smaller than for the Wrec distribntions (compare with the left fignre). 
The mass resolntions, derived from the Ganssian fits to the Wcorr distribntions in the 
region from /i —l.Scr to /i-|-1.3cr, are equal to 3.2 and 4.9 GeV, without and with overlaying 
events, respectively. Also the mean valnes obtained from the fit are closer to Mh than the 
valnes when MVec was nsed. 

The final Wcorr distributions for the signal and backgronnd events (with overlaying 
events included) are shown in Fig. 15.Ill For Mh = 120 GeV the most precise measurement 
of the Higgs-boson production cross section is obtained for the mass window between 108 
and 133 GeV, as indicated by arrows. In the selected Wcorr region one expects, after 
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one year of the Photon Collider running at nominal luminosity, about 4900 reconstructed 
signal events and 5400 background events {i.e. iisli^B ~ 0.9). This corresponds to the 
statistical precision of: 

A [r(;i o-7)br(a ^ bb)] 

r(ft ^ TT)BR(ft ^ 

The statistical precision calculated with the formula EZl should be considered a conser¬ 
vative estimate as it does not take into account our knowledge of the shape of the signal 
and background contributions. To exploit this additional information we can determine 
the measured number of signal events by the maximum likelihood method. In this pro¬ 
cedure the background contribution is assumed to be fixed and the likelihood function 
depends only on the total number of signal events. The Poissonian distribution of events 
in each bin is assumed. For the production of the Higgs boson with = 120 GeV 
we find that precision of the signal cross-section determination from the measured Wcorr 
distribution (without mass window cut) is 2.0%. 

The maximum likelihood method can also be used to estimate the systematic uncer¬ 
tainty of the cross section measurement. With this approach the systematic uncertainties 
of the estimated total background contribution and of the luminosity determination are 
included. We assume that both uncertainties are described by Gaussian distributions. 
The systematic uncertainty of the total background contribution is estimated to be about 
2%. This is based on the assumption that the background can be constrained from the 
dedicated run of the Photon Collider at lower value of ^/s^. If the center-of-mass-system 
energy is reduced by about 10 GeV, then only background events are measured (the Higgs 
boson will not be produced due to very low luminosity at = M^) while the detector 
performance remains almost unchanged. During half a year of the Photon Collider run¬ 
ning about 3000 events can be selected in the invariant-mass range corresponding to the 
optimal mass window for the Higgs-boson production measurement. The statistical un¬ 
certainty of 2% on the background contribution can be obtained which could be reduced 
by in a longer run if required. This uncertainty will result in the corresponding systematic 
error in the (independent) Higgs-boson production measurement. The J 2 = 0 luminosity 
contribution will be measured with precision of around 1% ua and this is assumed to be 
the expected uncertainty for the total luminosity (the \ Jz\ =2 luminosity contribution will 
be known with much better precision but it is small in the Higgs-resonance region). Using 
maximal likelihood procedure with assumed systematic uncertainties we obtain precision 
of 2.7% for < 7(77 h ^ hh) measurement at Mh =120 GeV. Therefore we can conclude 
that the systematic error of the measurement is of the order of 1.8%. 

The final result for r(/r — 77)BR(fi — bb) should be extracted from the measured 
event rate by applying correction for the selection efficiency. In our analysis the total 
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efficiency for signal events is only about 30%. The significant reduction of the signal 
events is due to the 6 -tagging cut. Therefore a very precise determination of flavour 
tagging efficiency will be crucial. To minimize influence of the uncertainties resulting 
from the Monte Carlo description of the detector performance we propose to use hadronic 
decays of Z bosons for 6 -tagging studies. After one year of the Photon Collider running at 
= 419 GeV, the sample of about 26000 77 —> ZZ events will be collected Taking 
into account branching ratios and the selection efficiency of about 80%, the expected 
number of Z ^ 66 decays will be about 5000. This sample will allow us to determine the 
flavour tagging efficiency for 6 -jets with statistical precision of about 1.4%. 

There is also an alternative solution: taking e+e“ data at = Mz for two to three 
months (the so-called GigaZ project), using the same detector and selection procedure 
as described above. Having about 10® events of Z-boson production we could determine 
flavour-tagging efficiencies with exceptional statistical precision of about 10“^. Also other 
systematic uncertainties could be significantly reduced with such a large sample of Z 
decays. Therefore we would like to stress that the possibility of operation in the 
collision point designed for the 77 mode should be guaranteed to make full use of the 
physics potential of the Photon Collider. 

For all methods, results for systematic uncertainties will require some extrapolation. 
However, we estimate that after a year of additional running, the total systematic uncer¬ 
tainty of the r(h — 77 )BR( 6 , — 66 ) measurement can be reduced to about 2 %. 

The background contribution due to 77 — W~^W~ production is included in the 
analysis for = 150 and 160 GeV. Hadronic decays of pairs result in 4-jet final 

state. However, in significant fraction of events only 3 or 2 jets are reconstructed. As 
cross section for production is very high, additional cuts, dedicated to suppress 

77 —background are introduced: 

1. Events are rejected if where is the highest invariant mass of 

the jet in the event. 

2. Total energy measured below 9tc^ Etc, should be less than Ce^c- 

3. Each jet in the event should contain at least tracks, i.e. > Cn^ where 

Aj™™ is the minimal number of tracks in jet for a given event. 

The first cut removes the events in which two jets from decay of one W boson have 
been joined into one jet by the jet-clustering algorithm. The second condition rejects 
events with substantial energy in the forward region, below 6tc, where one or two jets 
from W~^W~ decay could deposit their energy. The last cut, with the parameter value 
Cnt = 4, suppresses leptonic decays of W. Eor Mh = 160 GeV the optimal values for the 
two other thresholds are = 70 GeV and Cej,c = 90 GeV. 
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Figure 5.12: Statistical precision of r(h —> 77 )BR(/i —> bb) measurement for the SM 
Higgs boson with mass 120-160 GeV. Final results of this analysis are compared with our 
earlier results, which did not take into account some of the background contributions, 
distribution of the interaction point or overlaying events. The lines are drawn to guide 
the eye. 


We have performed a full simulation of signal and background events also for Mh = 
130, 140, 150 and 160 GeV choosing optimal e“e“ beam energy for each Higgs-boson mass. 
Statistical precision of F(h ^ 77 )BR(h —> bb) measurement was estimated in each case. It 
is equal to 2.1%, 2.5%, 3.4% and 7.7%, respectively. These results, together with the result 
for Mh = 120 GeV described above, are presented in Fig. 15.121 For comparison, our earlier 
results obtained without overlaying events, without various background contributions or 
without distribution of interaction point are also shown. For Mh = 160 GeV, after the full 
optimization of the selection cuts, better precision is obtained than in earlier analyses, 
which did not take into account some background contributions. 
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Chapter 6 


Production of heavy neutral Higgs 
bosons in the MSSM 


The analysis of the MSSM Higgs-bosons production closely follows the SM study described 
in Chapter El Therefore only these parts of the analysis which differ from the SM case 
are discussed in detail below. 

The analysis was developed for set / of the MSSM parameters, as defined in Section 
HU with tan /? = 7. To simplify the description the procedure is first presented for M 4 = 
300 GeV, and later extended to Ma = 200, 250 and 350 GeV. In each case the beam energy 
was chosen in such a way, as to provide the highest luminosity for 77 collisions with Jz = 0 
at = Ma- For detailed information about beam energies and 77 -luminosities, for the 
considered Higgs-boson masses, see Tab. 13.3.1 1 

For ~ 400 GeV about two 77 — hadrons events are expected on average in each 
bunch crossing. To suppress the influence of these events optimization of the 9tc cut was 
repeated, as described in AnnendixICland the value 9tc = 0.85, used in the SM analysis, 
turned out to be optimal also for the MSSM case. 

Selection criteria were optimized for the measurement of the higgs production cross 
section. The event distributions for | cos and \Pz\/E are shown in Fig. lfi.1l a,nd lfi.2[ 
respectively. For M 4 = 300 GeV the optimized cut values are Cq = 0.65 and Cp^ = 0.06 
(indicated by arrows). The measurement precision is estimated to be around 32% and 
22 %, after the | cos 6 *F*| cut and after the \Pz\/E cut, respectively. 

As in the SM case, bb events were selected by considering 6 -tag values for two jets 
with highest transverse momentum. For Ma = 300 GeV the selection region in the 
btag(jeti) ® btag(jet 2 ) plane which gives the best precision of the < 7(77 A, H ^ bb) 
measurement is shown in Fig. 16.31 Optimal bb selection, when overlaying events are taken 
into account, corresponds to the efficiencies Sh = 53%, Ebb = 47%, Ecc = 2.9%, and 
£uds = 0.5%, i.e. cc background suppression by a factor of Ebb/Scc ~ 16. This should be 
compared to results obtained when overlaying events are not included - then the optimized 
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Figure 6.1: Distributions of | cosOjetl’^^'^ for signal and background events. Only 77 — 
bb{g) events are shown for the background. The signal measurement precision for events 
fulfilling the cut | < Ce = 0.65, /a, is around 32%. 
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Figure 6.2: Distributions of \ Pz\/E for signal and background events. Only 77 — 

bb{g) events are shown for the background. The signal measurement precision for events 
fulfilling the cut \Pz\/E < Cp^ = 0.06, ^cr/a, is around 22%. 
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Figure 6.3: Ratio of ■yy ^ A, H —>■ bb events to 77 ^ bb{g), 77 —»• cc{g) and 77 ^ qq, 
q = u,d, s, events distributions in the plane btag(jeti) 0 btag(jet 2 ). The region is indicated 
by stars which gives the best precision measurement for <7(77 A, H ^ bb). 

selection corresponds to the efficiencies Eh = 57%, Ebf, = 52%, Ecc = 1-8%, and E^ds = 0.1%, 
i.e. Ebb!Sec ~ 29. This comparison shows that for heavy Higgs bosons overlaying events 
increase relative probability of cc-mistagging by a factor of about 2. Even larger effect is 
observed for light qnark-pair prodnetion as we observe increase of E^ds by a factor of 5 
despite the tighter selection cnt. 

In Fig. 16.41 we compare the (IFrec — 1 ^ 77 ) distributions, for signal events 77 —> % —> bb, 
before and after taking into account the overlaying events. Mass resolution, derived from 
the Ganssian fit in the region from /i —l.Scr to /x+1.3ct, is abont 8 and 12 GeV, respectively. 
The overlaying events and cut 9tc suppressing their contribution significantly influence 
the reconstructed mass distribution and result in a deterioration of the mass resolution. 
However, no significant drop in the selection efficiency is observed (the nnmber of selected 
events is abont 220 in both cases), in contrary to the SM case. Corresponding distribntions 
of {Wcorr — HTy^y) are compared in Fig. 16.41 The mass resolution for yy ^ A ^ bb events, 
for Ma = 300 GeV, is abont 7 GeV withont and abont 10 GeV with overlaying events. As 
the mass resolntion is mnch bigger than the mass difference between heavy nentral Higgs 
bosons H and A, Mh — Ma ~ 1.5 GeV, the separation of these two contributions will 
be impossible for the considered parameter set. Even for tan/3 = 3 the mass difference 
{Mh — Ma ~ 6.8 GeV) will be too small to resolve A and H contribntions, see Fig. 16.51 

The influence of overlaying events on the Higgs-boson production measurement is also 
illustrated in Fig. 16.61 and o where the Wrec and Wcorr distributions, respectively, are 
shown for signal and backgronnd events (for this comparison only quark-pair production 
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Figure 6.4: Mass resolution of the Higgs-boson A , Wrec — ^^77 (left), and Wcorr — 
(right), for the selected events for Ma = 300 GeV, obtained without and with overlaying 
events (OE). Mean and dispersion values, /r and a, from the Gaussian fit between /r —l.Sa 
and fi + 1.3cr are indicated. 
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Figure 6.5: Corrected invariant mass, Wcorr, distributions for the selected 77 ^ 66 

and 77 —A —66 events {Ma = 300 GeV), with overlaying events, for tan/3 = 7 (left) 
and tan/3 = 3 (right), respectively. For tan/3 = 3 also the sum of both contributions is 
shown. 
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Figure 6.6: Distribution of the reconstructed invariant mass, bFrec, for bb events (M^ = 
300 GeV) selected before (left plot) and after (right plot) taking into account overlaying 
events. Contributions of the H and A signal and of the quark-pair backgrounds are 
shown separately. Arrows indicate the mass windows optimized for the measurement of 
cr(77 A, H ^ bb). 
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Figure 6.7: As in Figs. 16.61 for the corrected invariant-mass, Wcorr, distributions. The 
statistical precision is 7.6% without overlaying events (left plot) and 9.4% with overlaying 
events (right plot). 
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Figure 6 . 8 : Distributions of the maximal jet mass, expected for signal and 

background contributions (overlaying events included). 


is taken as the background). In both cases distributions obtained without and with over¬ 
laying events are compared. We conclude that overlaying events significantly deteriorate 
the cross section measurement by increasing the selection efficiency of background contri¬ 
butions. As in the SM case, better estimates for measurement error Acr (77 —> A, —> 
66 )/cr (77 —> A, if —> bb) are obtained when using the variable Wcorr- 

For all considered heavy Higgs-boson masses in the MSSM, additional cuts to suppress 
77 —> W~^W~ background were also applied, as described for the SM Higgs boson with 
Mh = 160 GeV. For Ma = 300 GeV the optimal threshold values are CMj^t = 65 GeV and 
Cej,c = 80 GeV. The corresponding distribntions of events are shown in Fig. 16.81 and iGOl 
The resnlting estimate of the cross-section measnrement precision is abont 14%. After 
the additional cut on the number of tracks in each jet {Cnt = 4), the final result for 
Ma = 300 GeV is shown in Fig. 16.101 For the optimized Wcorr window the cross-section 
measnrement precision of 11 % is obtained. 

After discovery or a ’hint’ of the resonant-like excess of events at LHC or ILC the Pho¬ 
ton Collider can be used to confirm the existence and to measure the cross section for pro¬ 
duction of the new state. In this case the beam energy will also be optimized for the 
production of the observed resonance. Under this assnmption, we estimate precision ex- 
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Figure 6.9: Distributions of the energy below 6tc measured in the calorimeters, Eqal, 
for signal and background contributions, with overlaying events included. 


pected for 77 —> A, if ^ 66 cross section measnrement in the considered MSSM scenario 
after one year of experimentation. In Fig. 16.Ill the precisions for all considered values of 
Ma are shown and compared with our previous results. Thanks to the cuts optimization 
better estimated precision for Ma = 200 GeV has been obtained; earlier resnlts were 
obtained with cuts optimized for Higgs-bosons production for Ma = 300 GeV. 

Presented results were also used to estimate precision of the cross-section measurement 
at the Photon Collider for other parameter sets described in Section 14.11 For all sets the 
total widths of the Higgs bosons A and H are comparable or smaller than for the set 
/ and their masses are nearly degenerated (the mass difference is smaller or comparable 
with detector resolution). Precision estimates for < 7(77 ^ A, H ^ 66 ) measurement 
for all considered parameter sets, assnming tan/? = 7 are compared in Fig. 16.1 2\ The 
most precise measurement is expected for parameter sets II and III — precision hardly 
depends on Ma and is about 10%. The worst measurement is expected for scenario JV, 
i.e. the one considered by the CMS collaboration. 

For all considered values of Ma the dependence of the measurement precision on tan f3 
was studied and the results are shown in Fig. 16.131 The precision weakly depends on tan/? 
if parameter sets II or III are considered. In case of parameter sets I or IV the precise 


66 

















a 

3 

Qi 

O 

IT) 

U 

Qi 

Dh 

{/! 

a 

a> 

O) 

o 

O) 

s 

s 


180 


160 


140 


120 


100 


Ag/g =11.0% 



H+A signal 
M^=300 GeV 
Parameter set I 
tgP=7 

Background: 

H bb(g) 

□ cc(g) 

H W^W 
uu,dd,ss 

I-1 + - 

T T 


Total = 808 fb 


-1 


250 275 300 325 350 375 400 


W [GeV] 

corr ■- -■ 


Figure 6.10: Distributions of the corrected invariant mass, Wcorr, for signal and all 

considered background contributions, with overlaying events included. The best precision 
of 11 % for 77 A, if —> bb cross section measurement is achieved in the Wcorr window 
between 285 and 325 GeV. 


measurement will not be possible for low tan/I values, tan/I < 5. 

As already mentioned in Chapter El considered range of MSSM parameters covers the 
region of so-called LHC wedge where the expected number of reconstructed A and H 
production events will be too small to claim 5a discovery. Also the heavy MSSM Higgs- 
boson production signal at the ILC may not be sufficient. In such a case the Photon 
Collider was considered as the only machine which will be able to confirm the existence 
of the heavy neutral higgses. For M 4 = 300 GeV, assuming parameter set / and tan j3 = 
7, the expected statistical significance 6 of the A + H signal measurement at the Photon 
Collider is 5 = Hs/^/Wb = 7-2, in the optimal mass window. The significance of the 
Higgs-bosons measurement, 5, for different parameter sets considered in this analysis is 
shown in Fig. 16.141 The bands widths indicate the level of possible statistical fluctuations 
of the actual measurement: _ 


= i + ti 

a/Ab V Fb 

In three from the four considered MSSM scenarios, assuming tan P = 7, 5a discovery of 
heavy, neutral higgses will be possible in the whole range of M^, if only the optimum beam 
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Figure 6.11: Precisions of < 7(77 A, H ^ bb) measurement for MSSM parameter set 
I, for Ma = 200-350 GeV and tan/? = 7. Final results of this analysis are compared with 
our earlier results, which did not take into account some of the background contributions, 
77 —> W~^W~, distribntion of primary vertex or overlaying events. The lines are drawn 
to gnide the eye. 

energy is chosen (some hint for Ma value is expected from LHC and ILC measurements). 

For all considered values of Ma we also studied the significance of signal measurement 
as a fnnction of tan/?, for tan/? = 3-20. Resnlts obtained for different parameter sets 
are compared in Fig. 16.1 51 The estimated lower limit of the discovery region of LHC 
experiments (as presented by CMS collaboration jHj) is indicated by arrows. For all 
parameter sets the expected statistics of signal events for M 4 = 200-350 GeV will be 
sufficient to cover most of the considered MSSM parameters space. We can conclude that 
for Ma > 300 GeV the Photon Collider should be able to discover Higgs bosons for much 
lower values of tan/? than experiments at the LHC. 


68 





Figure 6.12: Precisions of < 7(77 ^ A, H ^ bh) measurement are shown for MSSM 

parameter sets I-IV, for Ma = 200-350 GeV and tan (3 = 7. The lines are drawn to guide 
the eye. 
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Figure 6.13: Precisions of cr (77 ^ A, H ^ bb) measurement are shown for = 200, 
250, 300 and 350 GeV, for MSSM parameter sets I-IV with tan/3 = 3-20. 
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Figure 6.14: Significances of < 7(77 ^ A, H ^ bb) measurement for MSSM parameter 
sets I-IV, for Ma = 200-350 GeV and tan/I = 7. The band widths indicate the level of 
possible statistical fluctuations of the actual measurement. 
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Figure 6.15: Significances of <7(77 A,H ^ bb) measurement for MSSM parameter sets 
I-IV, for Ma = 200, 250, 300 and 350 GeV, and tan/3 = 3-20. The band widths indicate 
the level of possible statistical fiuctuations of the actual measurement. 
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Chapter 7 
Summary 


Search for the Higgs boson and precise measurement of its properties are among most im¬ 
portant tasks of all future colliders. One of the measurements crucial for understanding of 
the Higgs sector and for the verification of the particle physics models is the measurement 
of r(/r — 77 ). The Photon Collider, which has been proposed as an extension of the 
linear collider project, is considered the best place to do this measurement. 

In this thesis the first fully realistic estimates for the precision of 77 —higgs 
bb cross-section measurement at the Photon Collider are obtained. Assumed photon 
collider parameters correspond to the superconducting Linear Collider project TESLA. 
The analysis is based on the realistic 77 -luminosity spectrum simulation. Due to the high 
beam intensity, resulting in high 77 -luminosity per bunch crossing, the contribution of 
overlaying events 77 —hadrons turns out to be sizable and affects event reconstruction. 
Crossing angle between beams resulting in the significant broadening of the interaction 
region is also taken into account. These two factors significantly affect the performance of 
the 6 -tagging algorithm. It is shown that the contamination of 77 —hadrons overlaying 
events in the signal can be reduced by rejecting low-angle tracks and clusters in the event. 
Additional cuts are proposed to suppress contributions from other background sources 
which are taken into account. For the realistic estimate of the heavy quark production 
background which is considered to be the most important one, the NLO QCD results are 
used. Other processes are taken into account using the LO cross sections implemented 
in PythiA: r-pair production, light-quark production 77 — qq, 77 —hadrons events 
(without directx direct processes and without other events in the bunch crossing), and 
77 — W~^W~ (with modified cross section to take into account polarization). 

After optimizing selection cuts and applying correction for escaping neutrinos from D- 
and 5-meson decays the quantity T{h ^ 77 )BR(/r — bb) for the SM Higgs boson with 
mass around 120 GeV can be measured with the precision of about 2 % already after one 
year of the Photon Collider running. The systematic uncertainties of the measurement are 
estimated to be of the order of 2%. The statistical precision of the measurement decreases 
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up to 7.7% for the SM Higgs boson with mass Mh = 160 GeV. For higher masses of the 
SM Higgs boson other decay channels are expected to give better precision of r(h ^ 77 ) 
measnrement, see e.g. |49| . 

The measurement discussed in this paper can be used to derive the partial width T{h ^ 
77 ), taking BR(h —^ bb) value from precise measurement at the International Linear 
Collider. With 2% accnracy on T{h 77 )BR(h — bb), as obtained in this analysis, and 
assuming BR(h —> bb) will be measured to 1.5% |1S|, Higgs-boson partial width T{h ^ 77 ) 
can be extracted with accuracy of about 2.5%. With this precision the measurement will 
be sensitive to the deviations from the SM coming from loop contribntions of new heavy 
charged particles. For example, heavy charged higgs contribntion in the SM-like 2HDM 
is expected to change V{h 77 ) by 5-10% |12j. Using in addition the result from the 
Linear Collider for BR(h — 77 ) (32], one can also extract F/j with precision of 10%. 

The presented analysis has been extended to the MSSM case by considering signal of 
heavy MSSM Higgs bosons A and H production with decays into bb. We study masses 
Ma = 200, 250, 300 and 350 GeV, and for each mass we choose an optimal e~e~ beam 
energy and selection cnts, as in SM case. Following |23|, we stndy MSSM parameter range 
which corresponds to the so-called “LHC wedge”, where A and H are almost degenerate 
in mass. 

Onr analysis shows that, for Ma ~ 300 GeV, the cross section for the MSSM Higgs- 
bosons production < 7(77 A, H ^ bb) can be measured with a statistical precision of 
about 11% already after one year of Photon Collider running. For other considered values 
of Ma it turns ont to be lower - from 16% to 21 %. Althongh this resnlt is less optimistic 
than the earlier estimate pg, still the photon-photon collider gives opportunity of a 
precision measurement of < 7(77 —>■ A, H ^ bb), assuming that we know the mass of the 
Higgs boson(s). Unfortnnately, separation of A and H will not be possible as the mass 
resolntion is abont 15 GeV. 

A discovery of MSSM Higgs-bosons requires energy scanning or a run with a broad 
Inminosity spectrnm, perhaps followed by the rnn with a peaked one. We estimate the 
significance expected for the Higgs prodnction measnrement for fonr different parameter 
sets and for tan (5 = 3-20. The discovery reach of the Photon Collider is compared with 
the estimated reach of the CMS experiment. For the optimnm energy choice the photon- 
photon collisions allow for the discovery of the Higgs bosons even for tan f3 valnes lower 
than the expected reach of the LHC experiments. Thus, at least partially, the Photon 
Collider will cover the so-called “LHC wedge”. For low tan/I values the measurement 
conld probably profit from nse of additional channels: decays of A and H to charginos 
and neutralinos, and in case of H also decays to hh. 

Presented analysis takes into account all corrections which were available to date. 
However, some aspects of the analysis can probably be improved in the fntnre. The first 
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one is the cut on the number of selected jets. In the current analysis 2 or 3 jets are required. 
These are the ’natural’ jet multiplicities for the signal. However, when full NLO generators 
with ’matrix element - parton shower matching’ will become available for the heavy 
quark background, then this study should be repeated and the cuts depending on the 
number of jets should be considered. Secondly, the jet-clustering algorithm with smaller 
Hcut values could be used. Narrower jets (possibly including tracks and clusters rejected 
by 6 *Tc-cut) could be used in the 4-jet kinematical fit to discriminate 77 — W^W~ 
background events more effectively. It is also not excluded that further minimization 
of the 77 —hadrons overlaying-events contribution is possible. We have not found any 
effect due to overlaying events on the spacial distributions of jet primary vertices, probably 
because tracks with largest impact parameters are rejected by the vertex fit. However, 
the cut on the impact parameter of single tracks and rejection of accompanying clusters 
could still be considered. 
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Appendix A 


Event generation with overlaying events 


In this Appendix we discnss the statistical principles of event generation, with emphasis 
on the possibility of having multiple events in single bunch crossing. We assume that for 
a given invariant mass of colliding photons, lEy.y, we know how to calculate the total cross 
section for the process a, aaiW^^). In the general case, for the 77 -luminosity spectrum 
(integrated over time), given by dL{W ^^)/the expected total number of events, Na, 
is: 

If we take the luminosity spectrum corresponding to one bunch crossing only then we 
obtain an average number of events of the process a per one crossing, jia- In the realistic 
event generation one must take into account the possibility that two or more 77 collisions 
in channel a take place in one crossing. Therefore in each generated event, corresponding 
to one bunch crossing, ka physical events of process a should be included, where ka is the 
random number from Poisson distribution with mean Ha' 

In this analysis the exact algorithm based on the Poisson distribution is used to generate 
77 —> hadrons events for which p ~ 1 . 

For small values of /Xq, /ia -C 1, an approximation is valid in which one generates Na 
bunch crossings with one event of the process a in each crossing, i. e. one assumes that the 
fraction of bunch crossings with exactly one a event is Ha- Actually, this fraction is smaller 
and is equal to P{ka = 1, /ia) = Moreover, in this simple approach bunch crossings 

with more than one a event are neglected; they contribute to the fraction P{ka > 1 ,/ia) 
of crossings. Inaccuracy of this procedure is of the order of Therefore it can be safely 
used for processes with pa 1. In this work this approach is used for following processes: 
77 — higgs 66 , 77 — QQ{g), 77 —> qq, 77 — W~^W~, 77 —r+r“ (for the process 
with the highest cross section, 77 — QQ{g), ja ~ 10 “^). 
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In general one would like to consider many independent physical processes. Here, to 
simplify the description, we assume that only three different processes have to be taken 
into account: a, b and c, with /i^, f^b < 0.01 and pc ~ 1. In the exact (but very ineffective) 
generation algorithm one would have to include ka, kb and kc events of corresponding 
processes for each bunch crossing where ki are random numbers coming from Poisson 
distributions P{ki,fii). However, for small Ha and fib one can neglect bunch crossings in 
which events of both process a and process b are present; the fraction of such crossings 
is P{ka > 0,fia)P{kb > 0,fib) ~ fJ^af^b < 10“"^. Therefore only bunch crossings with a and 
c, or b and c processes can be considered. As already pointed out, only bunch crossings 
with ka = 1 and kb = 1 can be taken into account in generation. In this case c events, 
with kc given by the Poisson distribution, are called overlaying events. However, one has 
to include also events with c process only. The fraction of bunch crossings with /Cc > 0 
and /Cq = /cb = 0 is equal to P{kc > 0, fi^P{ka = 0, fia)P{kb = 0, fib). As values of fia and 
fib are small P{ka = 0,fia)P{kb = 0,fib) ~ 1 , and it is enough to generate a sample of 
P{kc > 0, fic)Nc/fic bunch crossings with random number kc {kc > 0) of c events in each 
crossing. In this analysis processes 77 — higgs bb, 77 — W~^W~ and 77 — 
(corresponding to processes a or b) with overlaying events 77 — hadrons (process c) are 
generated according to this approach.^ 

Slightly different procedure can be applied if process a is a subclass of overlaying 
event process. Let us first consider the generation of processes a with overlaying events 
c in more detail. The fraction of bunch crossings in which there are two a events is 
P{ka = 2, fia) ~ /Ua/2. Usiug simple approach for generation of a events, i.e. generating 
Na bunch crossings with ka = 1, and overlaying a random number of c events, kc, one 
overestimates number of crossings with one a event by about /i^ and - by neglecting 
a possibility of having ka > I - underestimates number of bunch crossings with two a 
events by about /ia/2. However, there is also an alternative version of the generation 
algorithm. Instead of generating kc overlaying events from class c, one can generate ka+c 
overlaying events from joined class a + c where ka+c is a random number from distribution 
P{ka+c, f^a + Then, one obtains almost correct fraction of bunch crossings with /Ca = 1 
(about fia—fia) and overestimates a fraction of crossings with ka = 2 (by about fii/^). This 
approach is applied in this work for processes 77 — QQ{g) and 77 — qq (corresponding 
to the process a) with overlaying events 77 — hadrons {a + c) because the first two 
processes are included in generation of 77 — hadrons events.^ 


^ When the 77 ^ hadrons process is considered as a separate background contribution (only overlaying 
events) the direct processes 77 —> qq with q = u,d, s, c, b are excluded from generation as they are taken 
into account individually (77 ^ QQ{g) for Q = b,c and 77 ^ qq for q = u,d,s). 

^The processes 77 ^ QQig) are included in 77 ^ hadrons events generated with Pythia in the LO 
approximation for the process 77 ^ qq with q = c,b, and assuming unpolarized photon-photon collisions. 
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Appendix B 

Influence of higher order corrections 


In this Appendix we try to estimate the uncertainty of the presented results due to higher 
order corrections to the NLO QCD background process 77 — QQ{g)- As explained in 
Section P~2l the parton shower algorithm was not applied in the generation of heavy-quark 
background events as it would lead to double-counting of some gluon emissions. However, 
we can use it to estimate influence of higher order QCD corrections which are not included 
in the NLO generating program. The parton shower algorithm was applied to the QQ{g) 
final state generated with the NLO code, and the full analysis including cut optimization 
was repeated for the SM Higgs boson with mass of 120 GeV = 210 GeV) and for 

the MSSM Higgs bosons H and A for Ma = 300 GeV {^/^ = 419 GeV). However, here 
background contributions due to 77 —> r+r“ and resolved 77 —hadrons events were not 
included. The cut suppressing overlaying event contribution was not varied: cos 9 tc = 
0.85. As some gluon emissions are neglected in the generation without parton shower 
(our standard analysis), and with parton shower some emissions are double counted, we 
can expect that the actual event distributions will be somewhere in between. Thus, the 
difference between results obtained without and with parton shower can indicate the size 
of effects which are expected after higher QCD corrections are taken into account. 

In Table H results obtained with and without parton shower are compared. The 
precisions of cross section measurement for both cases are nearly the same. After applying 
the parton shower algorithm the 77 ^ bb{g) events have on average slightly lower 6 -tag 
values as i?-mesons are slower due to the additional gluon emissions. It can be clearly seen 
in the analysis of SM Higgs-boson production when in both cases we obtain similar optimal 
region of acceptance in the 6 -tag(jeti)G) 6 -tag(jet 2 ) plane (nearly the same efficiencies for 
higgs signal) but Ebb is smaller if the parton shower algorithm is applied. Moreover, the 
increase of cc-mistagging probability is observed. It is due to the fact that additional soft 
gluons broaden a c-jet and some tracks appear which are not present in events without 
parton shower. These modifications can cause the c-jet being sometimes similar to the 
6 -jet. 
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Signal process 

Aa/a [%] 

[%] 

£bb [%] 

£cc [%] 

Suds [%] 

210 

77 —h —66 

1.96 (2.03) 

57 (58) 

43 (50) 


0.16 (0.16) 

419 

yy ^ A, H ^ bb 

10.6 ( 11 . 0 ) 

59 (53) 

50 (47) 


0.6 (0.5) 


Table B.l: Results for precision measurement of signal jj ^ h ^ bb {Mh = 120 GeV) 
and ■yy ^ A,H ^ bb {Ma = 300 GeV), and optimal tagging efficiencies Eh, Ebb, £cc, and 
Euds- Values obtained with parton shower algorithm applied to QQ{g) parton system. 
In parentheses values for generation without parton shower are shown (our standard 
analysis). All cases were subject to the same procedure with cuts and higgs-tagging 
optimization. 


In case of the MSSM analysis the situation is more complicated as the wider optimal 
acceptance region in the 6 -tag(jeti)G)&-tag(jet 2 ) plane is chosen if parton shower is included 
(greater Eh). Then it is advantageous to take also signal events with jets having lower 6 -tag 
value because the 77 — bb{g) background events have slightly lower 66 -tagging efficiency. 

We conclude that our results are not sensitive to the higher order corrections for the 
heavy-quark background process under assumption than these corrections do not change 
significantly the overall normalization of the 77 — QQ{g) cross section. 
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Appendix C 

Studies on the influence of overlaying 
events 


In this Appendix the contribution of the 77 —> hadrons overlaying events to the hard 
77 scattering events in the detector is studied in detail. To understand better properties 
of overlaying events we first study particle distributions on the generator level. The 
expected distributions of the particle, energy and transverse energy flows in cos{9particle), 
where 6parUcie is the particle polar angle, are shown in Fig. IC.1 1 for hadron-likex hadron¬ 
like events at = 210 GeV. The average number of 77 —> hadrons events per bunch 
crossing, /r, is about 0.9. All distributions are strongly forward-peaked. This behaviour, 
expected for soft hadron-hadron interactions, is intensified in majority of events by the 
strong Lorentz boost of the 77 center of mass system. As shown in Fig. I(k2l and 1C.31 the 
same behaviour is observed for hadron-likex direct and directx direct processes for which 
/r ~ 0.1 and 0.01, respectively. Therefore, the cut on the particle polar angle seems to be 
the simplest way to reduce the contamination of particles coming from overlaying events. 
We studied the average contribution of overlaying events to the energy measured in the 
detector after a cut 9min on the minimal polar angle of particle. In Fig. 1C.41 and 1C.51 the 
probability distributions for the energy and trasverse energy flows are shown for two 9min 
values. The significant decrease of the average energy and transverse energy is clearly 
observed: the average energy changes from around 60 GeV without 9min cut to 20 and 
10 GeV for 9min = 80 and 250 mrad, respectively. 

After these promising results on the generator level, we studied the detector response 
to the overlaying events. The influence of low-angle particles can be reduced by ignor¬ 
ing tracks and clusters with polar angle less than 9tc- Following values of 9tc were 
considered: 131, 315, 450, 555, 643 and 722 mrad, corresponding to cos 9 tc values of 
(approximately) 0.99, 0.95, 0.90, 0.85, 0.80 and 0.75. The lowest value 9tc = 131 mrad 
corresponds to the expected size of the mask, 9mask- Probability distributions for re¬ 
constructed energy, Erec, for various values of 9tc are shown in Fig. 1C.61 and ICTlfor 
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hadron-like x hadron-like interactions 


hadron-like x hadron-like interactions 


hadron-like x hadron-like interactions 





Figure C.l: Distributions of particle, energy and transverse energy flow in cos{9particle) 
for hadron-likex hadron-like case for = 210 GeV. 
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Figure C.2: Distributions of particle, energy and transverse energy flow in cos{9parUcie) 
for hadron-likex direct case for ^/s^ = 210 GeV. 
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Figure C.3: Distributions of particle, energy and transverse energy flow in cos{9particle) 
for directx direct case for = 210 GeV. 
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Figure C.4: Energy distributions for 77 —hadrons event with 9min = 80 and 250 mrad. 


71 —> hadrons (|i=1.0 events per bunch crossing) yy ^ hadrons (p=1.0 events per bunch crossing) 




Figure C.5: Transverse energy distributions for 77 —hadrons event with 6min = 80 and 
250 mrad. 
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= 210 and 419 GeV, respectively. Probability distributions of all buch crossings 
(multiple overlaying events possible) are compared with probability distributions of buch 
crossing with only one 77 ^ hadrons event (visible in the detector). The low energy 
parts of the distributions are presented in detail in Fig. EJandEHl showing a dip due to 
the limited detector sensitivity around Erec ~ 0.5 GeV. The most significant reduction of 
77 —> hadrons contamination is obtained when increasing 9tc cut from 9mask to cos9tc = 
0.95 and 0.9. However, it is important to notice that the probability of observing no en¬ 
ergy in the detector {Erec < 0.2 GeV) does not increase significantly with increase of 9tc- 
Most events remain visible in the detector although with highly reduced energy. 

As the Pythia generator does not produce 77 —> hadrons events for < 4 GeV, we 
try to estimate the contribution of these events. In Fig. IG.lOl average energy reconstructed 
per one 77 — hadrons event, <Erec>, is shown as a function of for ^/s^e = 210 GeV 
and 419 GeV. Again, the strong reduction of reconstructed energy is observed with increase 
of 6 *rc- It is worth noticing that for =210 GeV the average reconstructed energy is 
slightly higher at given than for ^/s^e = 419 GeV. This is due to the larger boosts 
of the low energy 77 events at higher e“e“-beam energies. However, the reconstructed 
energy per bunch crossing is higher for = 419 GeV as the expected number of 

77 —> hadrons events increases from about 1 to about 2 per bunch crossing. From the 
average numbers of overlaying events per bunch crossing and the fact that the fraction 
of luminosity below = 4 GeV is equal approximately to 30%, we estimated that 
the neglected \ow-W^^ overlaying events would contribute to no more than 0.3 GeV (on 
average) to the energy measured in the detector if cos 9 tc ~ 0-85 is chosen. 

Particles coming from overlaying events, entering the detector with the angle larger 
than 9tc, are expected to form jets with small relative transverse momentum, /Et, 
where is the transverse momentum of the jet and E^ is the transverse energy of 
all energy-flow objects above 9tc- We studied reconstructed jet distributions for the 

^ h ^ bb and 77 —> A —events with overlaying events included, for various 
values of cos^T^c'-cut. For comparison the ’ideal’ measurement, i.e. without overlaying 
events (therefore using 9tc = 9mask) was also considered. For consistency with the final 
analysis we required that for the jet with the highest transverse momentum an angle 
cut was fulfilled: \cos9jet\ < 0.7, where 9jet is the polar angle of the jet. The pH^^/E t 
distributions for jets reconstructed in 77 —fi/A —> 66 events are shown for various values 
of cos9tc in Fig. IG.l ll ( h production for y/s77 = 210 GeV and Mh = 120 GeV) and in 
Fig. inHSl (A production for ^/s^ = 419 GeV and Ma = 300 GeV). For bunch crossing 
with 9tc = 9mask the contribution from overlaying events is clearly seen for p!^^/E t < 0 . 1 . 
Already for the cut cos 9 tc = 0.95 this contamination is significantly reduced and the 
distribution approaches the ’ideal’ one. 

With increasing the 9tc cut the influence of overlaying events on the Higgs-boson 
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Figure C. 6 : Probability distribution for reconstructed energy, -Free, for various values of 
9tc and for ^/s^ = 210 GeV. On the right plot the probability for a bunch crossing with 
only one 77 —> hadrons event which on the generator level has some particles above the 
detector mask. For clarity only centers of bins are connected by lines. 




Figure C.7: Probability distribution for reconstructed energy, E^ec, for various values of 
9tc and for = 419 GeV. On the right plot the probability for a bunch crossing with 
only one 77 —> hadrons event which on the generator level has some particles above the 
detector mask. For clarity only centers of bins are connected by lines. 
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production events is decreasing. Also the lower systematic uncertainty due to 77 — 
hadrons contributiuon coming from relatively large uncertainties of the cross sections for 
hadron-like interactions of photons can be expected. Hovewer, some tracks and clnsters 
coming from signal event can be ignored at the same time. This will decrease mass 
resolution and worsen precision of the measurement. In Fig. 1C.Ill and 1C.O distributions 
of the reconstrncted jet mass, Mj^t, are shown for different cos 9 tc cnts. We see that for 
cos 9Tc = 0.75 and 0.8 some jets loose most of their content and their mass is significantly 
underestimated {Mjet ^ 3 GeV). Distribution obtained for cos 9tc = 0.85 nearly recovers 
the ’ideal’ jet-mass distribntion having almost the same mean valne. 

Finally, we checked the infinence of the ^^c'-cnt on the Higgs-boson mass resolution. 
For this comparison the mass resolution was defined as the dispersion, a, of the Gaussian 
distribution fitted to IVrec — ^'y'y or Wcorr — W^j distribution. The fit was performed in the 
range (/r — a, ju + 1.5a), where fi is the mean of the Ganssian distribntion. As shown in 
Fig. lG.13l a.nd IG.lH mass resolutions do not change significantly for cos 9 tc < 0.9. Both, 
the obtained resolution values and the shape of distributions show that the cut values 
cos9tc = 0.99 and 0.95 are too weak to suppress the 77 —hadrons contribution. 

The full analysis (see GhapterEj) considering Higgs-boson production signal and the 
77 —> QQ{g) background was performed for processes 77 —h —66 (for Mh = 120 GeV) 
and 'j'j ^ A, H ^ bb (for Af4 = 300 GeV), for all considered values of cos9tc- The 
best precision of the cross-section measnrement was obtained for cos 9tc = 0.85 and 
0.9. Because higher 9tc cut does not deteriorate the precision, but has the advantage of 
redncing the systematical nncertainty dne to 77 — hadrons events, the valne cos 9 tc = 
0.85 was chosen as the optimal one. 
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Figure C. 8 : Probability distribution for reconstructed energy, -Free, for various values of 
9tc and for ^/s^ = 210 GeV. On the right plot the probability for a bunch crossing with 
only one 77 —> hadrons event which on the generator level has some particles above the 
detector mask. For clarity only centers of bins are connected by lines. 




Figure C.9: Probability distribution for reconstructed energy, E^ec, for various values of 
9tc and for = 419 GeV. On the right plot the probability for a bunch crossing with 
only one 77 —> hadrons event which on the generator level has some particles above the 
detector mask. For clarity only centers of bins are connected by lines. 
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Figure C.IO: Average reconstructed energy, <Erec>, dependence on 77 invariant mass, 
for bunch crossings with one 77 —hadrons event. Various values of 9tc for ^/^ = 
210 GeV (left) and 419 GeV (right) are considered. 
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Figure C.ll: Distribution of /Et (left) and Mj^t (right) for jets reconstructed in 
■yy ^ h ^ bb event, for ^/s^ =210 GeV and Mh = 120 GeV. Various values of 9tc are 
considered. For comparison results for the ideal case without overlaying events are also 
shown (no OF). The mean value of the Mjet distribution is presented for each case. 




Figure G.12: Distribution of /Et (left) and Mj^t (right) for jets reconstructed in 
yy ^ A ^ bb event, for = 419 GeV and Ma = 300 GeV. Various values of 9tc are 
considered. Also results for the ideal case without overlaying events are shown (no OE). 
The mean value of the Mjet distribution is presented for each case. 
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Figure C.13: Distribution of Wrec — (left) and Wcorr — fFy-y (right) for accepted 
77 —> h —> 66 events, for ^/s^ = 210 GeV and = 120 GeV. Various values of 9tc are 
considered. Also results for the ideal case without overlaying events are shown (no OE). 
The mean valne, p, and the dispersion, a, of the fitted Ganssian distribntion is presented 
for each case. The fit was performed in the range (fu. — a, /i + l.Sa). 




Figure C.14: Distribution of Wrec — W^^ (left) and Wcorr — 1^77 (right) for accepted 
■yy ^ A ^ bb events, for ^/s^ = 419 GeV and M 4 = 300 GeV. Various values of 9tc are 
considered. Also resnlts for the ideal case withont overlaying events are shown (no OE). 
The mean valne, /i, and the dispersion, a, of the fitted Ganssian distribntion is presented 
for each case. The fit was performed in the range (/i — a, /i + l.Sa). 
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Appendix D 
Orlop 


In this Appendix we describe the package Orlop (OveRLaying events fOr Photon collider) 
which was developed to take 77 —> hadrons overlaying events into account in our analysis. 
The package includes subroutines which read in the luminosity spectra for 77 collisions, 
appropriately rescale energies of beam photons (if necessary), generate 77 —> hadrons 
events with variable beam energies, and add random number of the overlaying events to 
the earlier generated events of the studied hard scattering process. 

Luminosity spectra obtained from the full simulation of 77 collisions PI are used to 
account properly for the low energy part of the spectrum. The results of the simulations 
are available only for three values of e“e“ invariant mass, ^/s^ = 200, 500 and 800 GeV. 
The luminosity spectrum for user-dehned value of a/s^ is obtained by rescaling photon 
energies from the closest set by the factor where and are 

maximal energies of photons coming from electrons scattering off one laser photon for 
a/s^ and respectively. To obtain and values the CompAZ program is 

used The Orlop package can be easily extended to provide overlaying events also 
to collisions if the beamstrahlung photon spectrum is implemented. 

The main program gen_orlop, running as a separate process, uses Pythia to generate 
77 —> hadrons events with variable beam energies, according to the realistic luminosity 
spectrum. As Pythia sometimes does not succeed to hadronize events with very low in¬ 
variant mass, gen_orlop uses only 77 collisions with = 4 GeV (see Appendix 

o for detailed discussion). The program sets following values for Pythia steering vari¬ 
ables: MSTP(14)=10 to properly include all kinds of real photon interactions, and MSEL=2 
to switch on elastic and diffractive processes. The crossing angle and the interaction point 
distribution are taken into account. At the beginning the program is run for 50000 events 
to calculate more precisely internal Pythia weights which influence contributions of var¬ 
ious event classes. The calculated total cross section is used to determine the average 
number of overlaying events per bunch crossing, ja. 

The subroutine add_orlop is called to include overlaying events during generation 
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of the considered processes. This subroutine adds particles which are provided by the 
program gen_orlop to the Pythia event record. The call sequence of add_orlop is 
shown in the following listing where only parts relevant for use of Orlop are present. 

* Example of ORLOP use: 

* Higgs-boson production ’gamma gamma -> h -> b bbar’ 

* with overlaying events ’gamma gamma -> hadrons’, 
program produce_higgs_with_oe 

* Set Higgs-boson production ’gamma gamma -> h -> b bbar’ 

* Initialization of ORLOP 

call add_orlop(l,ibeams,sqrts_ee,n_oe_per_bc) 

* Additional information (optional call) 

call add_orlop(2,n_collisions,luminosity,cross.sec) 

* Loop over bunch crossings 
do ibc=l,n_bc 

* Set photons energies for ’gamma gamma -> h -> b bbar’ 

* Generate ’gamma gamma -> h -> b bbar’ event 
call pyevnt 

* Add overlaying events ’gamma gamma -> hadrons’ 
call add_orlop(3,n_oe,e_oe,pz_oe) 

enddo 

* End of loop over bunch crossings 

* Stop ORLOP 

call add_orlop(4,n_oe_tot,n_rewind,cross.sec) 
end 


In calls of the subroutine add_orlop(if lag, ivar ,varl ,var2) the action undertaken 
and also the meaning of variables ivar ,varl ,var2 depend on the value of the first pa¬ 
rameter if lag as explained below. 

The package Orlop, which contains also a manual, shell scripts for convenient run¬ 
ning of programs, and two examples of generator-level analysis, will be made publicly 
available after publication of results obtained in this thesis. 
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iflag=l 


Initialization. The main parameters for 77 —hadrons generation are passed 
to the gen_orlop program. Variables ivar and varl shonld be eqnal to the 
beam type and respectively. These - together with if lag - are the 

only input parameters. Beam type is at the moment defined only for 77 
collisions in the Photon Collider at TESLA. Variable var2 is set to the 
average number of overlaying events per bunch crossing, /i, calculated by 
gen_orlop. 

if lag=2 Optional call allowing to obtain additional information from Orlop. On 

return, variable ivar is equal to the number of beam collisions which are 
available in the simulated set. Variables varl and var2 are equal to the 
luminosity of used spectrum [pb“^s“^] and to the 77 —hadrons cross section 
[pb]. 

iflag=3 The call which actually adds overlaying events to the event record. The 
returned value of variable ivar is equal to the number of 77 —hadrons 
events added. This is a random number from Poisson distribution with 
mean jj, (see description of if lag=l). Variables varl and var2 are set to the 
total energy and to the longitudinal momentum of overlaying events. Before 
overlaying events are added all undetectable objects (except neutrinos) are 
removed from event record with call pyedit(l). 

if lag=4 Terminate 77 —hadrons generation by stopping gen_orlop. Summary of 

overlaying events generation is returned. Variable ivar is set to the number 
of all added 77 —hadrons events. Variable varl is equal to the number of 
times the available set of beam collisions was used. Variable var2 is equal to 
the 77 —hadrons final value of the cross section [pb] calculated by Pythia; 
this can be compared to the value provided with if lag=2 to check stability 
of the result. 
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